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Abstract: Using holographic entanglement entropy for strip geometry, we construct 
a candidate for a c-function in arbitrary dimensions. For holographic theories dual to 
J> Einstein gravity, this c-function is shown to decrease monotonically along RG flows. 

A sufficient condition required for this monotonic flow is that the stress tensor of the 
O matter fields driving the holographic RG fiow must satisfy the null energy condition 

^ over the holographic surface used to calculate the entanglement entropy. In the case 

where the bulk theory is described by Gauss-Bonnet gravity, the latter condition alone 
^ is not sufficient to establish the monotonic flow of the c-function. We also observe 

that for certain holographic RG flows, the entanglement entropy undergoes a 'phase 
transition' as the size of the system grows and as a result, evolution of the c-function 
may exhibit a discontinuous drop. 
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1. Introduction 

Zamolodchikov [1] showed that renormalization group (RG) flows of two-dimensional 
quantum field theories were governed by a remarkable underlying structure. One im- 
portant feature was that there exists a positive definite function C2, which decreases 
monotonically along the RG flows. At the fixed points of the RG flow, this function is 
stationary and coincides with the central charge c of the conformal field theory (CFT) 
describing the fixed point. A direct consequence for any RG flow connecting two such 
fixed points is then that 

[^].v > ['L- (1-1) 
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More recently, Casini and Huerta [2] developed an elegant reformulation of Zamolod- 
chikov's c-theorem in terms of entanglement entropy in two dimensions. In their con- 
struction, the c-function was defined as 

where S{i) denotes the entanglement entropy for an interval of length i. Then it follows 
that dc2/di < from the strong subadditivity property of entanglement entropy, as 
well as the Lorentz symmetry and unitarity of the underlying quantum field theory 
(QFT). Therefore, as the QFT is probed at longer distance scales, i.e., one increases 
(., this c-function (1.2) decreases monotonically. Further, for a two-dimensional CFT, 
the entanglement entropy is given by [3, 4] 

^cPT = ^log(£/(5) + c', (1.3) 

where c is the central charge, (5 is a short-distance regulator and c' is a non-universal 
constant (independent of i). Hence C2 = c at RG fixed points. 

As a generalization of the two-dimensional c-theorem, Cardy [5] conjectured that 
the central charge associated with A-type trace anomaly - see eq. (2.8) - should decrease 
monotonically along RG flows for QFT's in any even number of dimensions. Of course, 
in two dimensions, this proposal coincides precisely with Zamolodchikov's result (1.1) 
since c = 12 A. Cardy's conjecture was extensively studied in d = A and a great deal 
of support was found with nontrivial examples, including perturbative fixed points [6] 
and supersymmetric gauge theories [7, 8, 9].^ Recently, a remarkable new proof of this 
c-theorem was presented for any four-dimensional RG flow connecting two conformal 
fixed points [12]. This result draws on earlier work involving the spontaneous breaking 
of conformal symmetry [13] and bounds on couplings in effective actions [14]. It remains 
to determine, however, how much more of the structure of two-dimensional RG flows 
carries over to higher dimensions.^ 

As we will review below, support for Cardy's generalized c-theorem was also estab- 
lished using the AdS/CFT correspondence [18, 19, 20]. One of the advantages of the 
investigating RG flows in a such holographic framework is that the results are readily 
extended to arbitrary dimensions . In particular then, the analysis of holographic RG 
flows identified a certain quantity satisfying an inequality analogous to eq. (1.1) for any 



^Note that a flaw was recently found [10] in a proposed counter-example [11] to Cardy's conjecture, 
related question which has seen active discussion in the recent literature is whether or not 
there exist interesting QFT's in higher dimensions which exhibit scale invariance but not conformal 
invariance [15, 16]. Of course, in two dimensions, it is proven that scale invariant QFT's are also 
conformally invariant [17]. 
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dimension, that is, for both odd and even numbers of spacetime dimensions. Since the 
trace anomaly is only nonvanishing for even d, a new interpretation was required for 
odd d. Ref. [20] identified the relevant quantity as the coefficient of a universal contri- 
bution to the entanglement entropy for a particular geometry in both odd and even d. 
These holographic results then motivated a generalized conjecture for a c-theorem for 
RG flows of odd- and even-dimensional QFT's. For even d, this new central charge was 
shown to precisely match the coefficient of the A-type trace anomaly [20] and so this 
conjecture coincides with Cardy's proposal. For odd d, it was shown that this effective 
charge could also be identified by evaluating the partition function on a rf-dimensional 
sphere [21] and so the conjecture is connected to the newly proposed F-theorem [22]. 

The above developments motivated the present paper which examines the the con- 
nections between entanglement entropy and RG flows in a holographic framework. 
Earlier work in this direction can be found in [23, 24, 25]. Here, we make a simple 
generalization of the c-function in eq. (1.2) to higher dimensions and then use a holo- 
graphic framework to examine its behaviour in RG flows. We are able to show that 
subject to speciflc conditions, the flow of the c-function is monotonic for boundary 
theories dual to Einstein gravity. In examining speciflc flow geometries, we also flnd 
that the entanglement entropy undergoes a 'flrst order phase transition' as the size of 
the entangling geometry passes through a critical value. That is, in our holographic 
calculation, there are competing saddle points and the dominant contribution shifts 
from one saddle point to another at the critical size. 

An overview of the paper is as follows: In section 2, we review the standard deriva- 
tion of holographic c-theorems with both Einstein gravity and Gauss-Bonnet gravity 
in the bulk. We stress that in either case, the monotonic flow of the c-function requires 
that the matter flelds driving the holographic RG flow must satisfy the null energy 
condition. In section 3, we discuss the holographic entanglement entropy for the 'strip' 
or 'slab' geometry and construct a c-function which naturally generalizes eq. (1.2) to 
higher dimensions. In section 4, we show that for an arbitrary RG flow solution in 
Einstein gravity, this c-functions decreases monotonically if the bulk matter flelds sat- 
isfy the null-energy condition. Section 5 considers explicit examples of holographic RG 
flows and demonstrates that in certain cases, the entanglement entropy undergoes a 
'phase transition.' As a result, the c-function exhibits a discontinuous drop along these 
RG flows. In section 6, we examine holographic RG flows with Gauss-Bonnet gravity 
and there, we flnd that the null-energy condition is insufficient to constrain the flow of 
our c-function to be monotonic. We conclude with a brief discussion of our results and 
future directions in section 7. Appendix A presents certain technical details related 
to the discussion in section 4. In the appendix B, we discuss holographic RG flow 
solutions in Gauss-Bonnet gravity. Finally, appendix C describes the construction of a 



- 3 - 



bulk theory for which the holographic flow geometries examined in section 5 would be 
solutions of the equations of motion. 

While we were in the final stages of preparing this paper, we learned of a similar 
study of entanglement entropy and holographic RG flows appearing in [26]. 



2. Review of holographic c-theorems 

Here we begin with a review of the holographic c-theorem as originally studied by 
[18, 19] for Einstein gravity. These references begin by constructing a holographic 
description of RG flows. The simplest case to consider is ((i+l)-dimensional Einstein 
gravity coupled to a scalar fleld: 



J / d'+'x ^(R- {d<f>y - (2.1 



21^-^ J V 2 

We assume that the potential V{(j)) has various critical points where the potential 
energy is negative, i.e., 



d{d-l) 2 , 
= - \2 ^here — 



0. (2.2) 



Here L is some convenient scale while the dimensionless parameters distinguish the 
different flxed points. At these points, the gravity vacuum is simply AdS^+i with the 
curvature scale given by = L"^ /af. 

Now in the context of the AdS / GET correspondence, the bulk scalar above is dual 
to some operator O and the flxed points (2.2) of the scalar potential represent the 
critical points of the boundary theory. In particular then, with an appropriate choice 
of the bulk potential, O will be a relevant operator for a certain flxed point and so 
an RG flow will be triggered by perturbing the corresponding critical theory by this 
operator in the UV. Of course, the holographic description of this RG flow is that the 
scalar fleld acquires an nontrivial radial proflle which connects two of the critical points 
in eq. (2.2). The bulk geometry for this solution can be described with a metric of the 
following form [18, 19] 

ds^ = e2^(") r/,j dx'dx^ + dr^ . (2.3) 

Here, the radial evolution of the geometry is entirely encoded in the conformal factor 
A{r). At a flxed point where the geometry is AdS^+i, the conformal factor is simply 
A{r) = r /L where again L is the AdS curvature scale. Implicitly, we will assume that 
asymptotic UV boundary is at r oo while the IR part of the solution corresponds 
to r — > — oo. Hence for an RG flow between two flxed points as described above, the 
metric (2.3) approaches that of AdS^+i in both of these limits. 
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Now following [18, 19], we define: 



o^dir) = ^ , (2.4) 

^ ' T{d/2) (i^A'ir))''-^ ^ ' 

where 'prime' denotes a derivative with respect to r. Then for general solutions of the 
form (2.3), one finds 



T{d/2) ii-^A'{r) 
T{d/2) ii-^A'{r) 



-dir) = - ^.Z:/.,,^,, A"ir) (2.5) 

(T\ - T\) > . 



Above in the second equality, Einstein's equations were used to eliminate A"{r) in 
favour of components of the stress tensor.^ The final inequality assumes that the matter 
fields obey the null energy condition [27]. Now given the usual connection between r 
and energy scale in the CFT, eq. (2.5) indicates that a(r) is always increasing as we 
move from low energies to higher energy scales. Further, if the flow function (2.4) is 
evaluated for an AdS background, one finds a constant: 

oZ = (idir)\ = ^, , , , , -I . (2.6) 
d dK yi^^s T{d/2) £f-^ ^ ' 

Hence if we compare this constant for the UV and IR fixed points of the holographic 
RG flow, we find the holographic c-theorem: 

hS]„v > Hl^ (2-7) 
To make closer contact with the dual CFT, we recall the trace anomaly [28, 29], 

{T\)=Y,BnIn-2 {-f/^A Ed , (2.8) 

n 

which defines the central charges for a CFT in an even number of spacetime dimensions. 
Each term on the right-hand side is a Weyl invariant constructed from the background 
geometry. In particular, is the Euler density in d dimensions while the J„ are 
naturally written in terms of the Weyl tensor (as well as its covariant derivatives), 
e.g., see [30]. Note that in eq. (2.8), we have ignored the possible appearance of a 
conformally invariant but also scheme-dependent total derivative. 



^Note that for the scalar field theory in eq. (2.1), we have T*t - T''^ = -(0')2/2 < 
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A holographic description of the trace anomaly was developed [31] and can be 
applied to the AdS^+i stationary points in the present case (for even d). These calcula- 
tions show that a^, the value of the flow function at the fixed points, precisely matches 
the A-type central charge in eq. (2.8), i.e., a*^ = A for even d [20]. Hence with the as- 
sumption that the matter fields obey the null energy condition, the holographic CFT's 
dual to Einstein gravity satisfy Cardy's conjecture of a c-theorem for quantum field 
theories in higher dimensions [5]. Of course, one must add the caveat that for these 
holographic CFT's, i.e., those dual to Einstein gravity, all of the central charges in 
eq. (2.8) are equal to one another [31]. Hence the holographic models (2.1) considered 
above can not distinguish between the behaviour of A and Bn in RG flows. 

It has long been known that to construct a holographic model where the various 
central charges are distinct from one another, the gravity action must include higher 
curvature interactions [32]. In part, this motivated the recent holographic studies of 
Gauss-Bonnet (GB) gravity [33] — for example, see [34]. In section 6, we will extend 
our discussion of holographic RG flows to GB gravity with the following action 



'-^J <^'^'-^9 {R+ ^,^ . 3) ^. ^ -^O^Y - V(^) j . (2^9) 
where 



— RabcdR"" — 4:RabR"^ +R ■ (2.10) 

As before, we again assume the scalar potential has various stationary points as in 
eq. (2.2), where the energy density is negative. Note that for convenience, we are 
using the same canonical scale L which appears for the critical points in eq. (2.2) in 
the coefficient of the curvature-squared interaction in eq. (2.9). Hence the strength of 
this GB term is controlled by the dimensionless coupling constant, A. We write the 
curvature scale L of the AdS vacuum as = L^//oo where the constant /oo satisfies 
[20] 

a- -/oo + A/^ = 0. (2.11) 
In general, eq. (2.11) has two solutions but we only consider the smallest positive root 



(2.12) 

with which, in the limit A — )■ 0, we recover f^o = and = L^/a], as discussed 
above for Einstein gravity. One would find that graviton fiuctuations about the AdS 
solution corresponding to the second root are ghosts [35, 36] and hence the boundary 
CFT would not be unitary. The theory (2.9) is further constrained by demanding that 
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the dual boundary theory respects micro-causahty or alternatively, that it does not 
produce negative energy fluxes [34, 37]. 

For our present purposes, the most important feature of GB gravity (2.9) is that the 
dual boundary theory will have two distinct central charges. To facilitate our discussion 
for arbitrary > 4, we would like to define two central charges that appear in any CFT 
for any d - including odd d - and hence for our pursposes, the trace anomaly is not a 
useful definition of the central charges. Following [37, 20], we consider: 

d-l 

CT = ^irru^\^\ [1-2A/J, (2.13) 





[- 


T{d/2) 








r(rf/2) 


Up 



d-l 



(2.14) 



The first charge Ct is that controlling the leading singularity of the two-point function 
of the stress tensor.^ The second central charge can be determined by calculating 
the entanglement entropy across a spherical entangling surface [20]. Using the results 
of [39], it was further shown [20] that is the central charge appearing in the A-type 
trace anomaly in even dimensions, i.e., a*^ = A in eq. (2.8). In terms of these central 
charges, the micro-causality constraints, referred to previously, are conveniently written 
as [37] 

<^<1 (2.15) 



d(d - 2) - 2 - a* "2 

Now assuming the existence of bulk solutions describing holographic RG flows for 
the GB theory (2.9),^ we can establish a holographic c-theorem following the analysis 
of [20]. We begin by constructing two flow functions [20]: 

^-"•'" rW2) (UA\r)f-^ (l-2ALM-(r)-). (2,16) 

"^<^' ^ fS) T^7W^^^ k - ^^^"'^■^'<^)^) ■ P-"' 

These expressions were chosen as the simplest extensions of eq. (2.4) which yield the 
two central charges above at the fixed points, i.e., ad{r)\Ads = 0^ and CT{T)\Ads = Ct 
— recall that Air) = r/L for the AdS vacua. Now let us examine the radial evolution 



^Here, as in [38], we have normalized Ct so that in the Hmit A — > 0, Ct = aj. This choice is 
shghtly different from that originally presented in [37], i.e., Ct|[37] = 331 ^-^^r^ CTlhere- 
^Appendix B includes a discussion of one approach to constructing such solutions. 
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of ad{r) in a holographic RG flow: 

a'(r) = id-l)7r''/^ A"{r) (l ~ 2XL' A'(rf] (2.18) 



T{d/2)i^-^A'{rY 



> 0. 



Here, the equations of motion for GB gravity (see eq. (B.12)) have been used to trade 
the expression in the first line for the components of the stress tensor appearing in 
the second line. As before with Einstein gravity, we assume the null energy condition 
applies for the matter fields for the final inequality to hold. In eq. (2.9), the matter 
contribution is still a conventional scalar field action and so just as before Tt — T/ = 
— (0')^/2 < 0. With this assumption, it then follows^ that ad{r) evolves monotonically 
along the holographic RG flows and we can conclude that the central charge is always 
larger at the UV fixed point than at the IR fixed point. Hence we recover precisely the 
same holographic c-theorem found previously with Einstein gravity, namely. 

One can also consider the behaviour of Ct along RG flows 

r{d/2)£t^A'{rY ^ ^ V d 

but there is no clear way to establish that Crir) has a definite sign. Hence this 
holographic model (2.9) seems to single out as the central charge which satisfies a c- 
theorem. This result has also been extended to holographic models with more complex 
gravitational theories in the bulk:^ quasi-topological gravity [20], general Lovelock 
theories [41, 42], higher curvature theories with cubic interactions constructed with the 
Weyl tensor [20] and f{R) gravity [41]. The result is also established for holographic 
models where the RG flow is induced by a double-trace deformation of the boundary 
CFT [43]. Given the relation ajj = A in even dimensions, these holographic results 
support Cardy's proposal [5] that the central charge A (rather than any other central 
charge) evolves monotonically along RG flows. However, it is even more interesting 
that these results suggest that a similar behaviour also occurs for the central charge 
in odd dimensions. Further while the original field theory definition of involved a 
calculation of entanglement entropy [20], it was shown that the same charge can also 



Cr'ir) = - ^M^;JZ.. A"{r) ( 1 - 2^ A A'irf) (2.20) 



^We note that some additional arguments are needed to ensure that there are no problems with 
A'{r) < for odd d [20]. 

^Similar results were also found to apply in the context of cosmological solutions [40]. 
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be identified by evaluating the partition function on S"^ [21]. Hence tlie exciting new 
field theoretic results of [22] provide further evidence for the same c-theorem in odd 
dimensions. 

In any event, a key requirement for the holographic c-theorem to hold for Einstein 
gravity (2.7) or for GB gravity (2.19) is that the matter fields obey the null energy 
condition. Of course, this holds when these gravitational theories are coupled to a sim- 
ple scalar field, as in eqs. (2.1) and (2.9), this constraint is trivially satisfied. However, 
phrasing the constraint in terms of the null energy condition allows for more general 
scenarios for the matter fields driving the holographic RG flow. We should add that the 
same constraint also ensures the holographic c-theorem holds for all of the extensions 
of the bulk gravity theory mentioned above. We might mention that violations of the 
null energy condition quite generally lead to instabilities [44] and so it is a natural con- 
straint to define a reasonable holographic model. In the following, we will also see that 
the same constraint can be related to the monotonic flow of a holographic c-function 
defined in terms of entanglement entropy. 



3. Holographic entanglement entropy and a c-function 

Before beginning our holographic analysis, we must first identify a candidate c-function 
using entanglement entropy for d > 3. Recall that [2] identifies such a c-function for 
two-dimensional quantum field theories as 

where S is the entanglement entropy for an interval of length i on an infinite line. As 
described above, using the result (1.3) for the entanglement entropy oid = 2 CFT's, one 
finds C2 = c at any fixed points of the RG flows, i.e., eq. (3.1) yields the central charge 
of the underlying CFT at the fixed points. We would like to emulate this construction 
in higher dimensions. However, one should recall that in general the entanglement 
entropy for field theories in higher dimensions will contain many (non-universal) power 
law divergences depending on the geometry of the entangling surface, e.g., see eq. (7.5). 
Hence we expect a simple derivative with respect to some scale characteristic of the 
entangling surface will typically yield a result which depends on the cut-off. While there 
may be various strategies to avoid this outcome - see further discussion in section 7 - 
here we take the following simple approach: First we note that, at the fixed points, the 
power law divergences are geometric in origin and all but the leading area-law terms 
vanish if the geometries of the background and the entangling surface are both flat. 
Hence we consider a 'strip' or 'slab' geometry, where the entangling surface consists of 
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two parallel flat ((i-2)-dimensional planes separated by a distance £ in a flat background 
spacetime. The entanglement entropy (of a CFT) then takes the simple form [45, 46] 

where ad and (3^ are dimensionless numerical factors and H is a(n infrared) regulator 
distance along the entangling surface - we assume that H ^ i. That is, H'^~^ is the area 
for each of the planes comprising the entangling surface and so the flrst contribution in 
eq. (3.2) is simply the usual area law term. The coefficient of the second finite term is 
proportional to a central charge in the underlying d-dimensional CFT, which we denote 
Cd- Hence we can isolate this central charge by writing 

£d-l 

Cd = ■ (3.3) 

Hence we are naturally lead to consider the quantity 

£d-l 

Cd = Pd ^ (3.4) 

as a candidate for a c-function along the RG fiows, so that q = Cd at the fixed points 
of the fiow. We will identify the precise value of the coefficient (3d with our holographic 
calculations below — see eq. (3.11). Comparing eqs. (3.1) and (3.4), we can view the 
latter expression as the simplest generalization of the two-dimensional c-function (3.1) 
to higher dimensions. At the outset, we wish to say that we will find below that will only 
be able to prove that this candidate c-function actually decreases monotonically along 
RG fiows for holographic models with Einstein gravity in the bulk. However, another 
goal in the following analysis is to connect the behaviour of this c-function defined 
using holographic entanglement entropy with the standard discussions of holographic 
c-theorems [18, 19, 20]. We should also mention that eq. (3.4) was previously suggested 
as a c-function in [46]. 

3.1 Holographic entanglement entropy on an interval 

In this section, we derive some of useful results to evaluate eq. (3.4) for holographic 
RG fiows in following sections. The holographic models in sections 4 and 5 will be 
described by Einstein gravity in the bulk, while we will consider GB gravity [33] in 
section 6. 

The seminal work of Ryu and Takayanagi [45, 46] provided a holographic construc- 
tion to calculate entanglement entropy. In the ci-dimensional boundary field theory. 
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the entanglement entropy between a spatial region V and its complement V is given 
by the following expression in the ((i+l)-dimensional bulk spacetime: 

SiV) = -^^ min[A{v)] (3.5) 

where v ^ V indicates that f is a bulk surface that is homologous to the boundary region 

V [47, 48]. In particular, the boundary of v matches the 'entangling surface' dV in the 
boundary geometry. The symbol 'min' indicates that one should extremize the area 
functional over all such surfaces v and evaluate it for the surface yielding the minimum 
area.^ Eq. (3.5) assumes that the bulk physics is described by (classical) Einstein 
gravity and we have adopted the convention: ip~^ = SttG^ Hence the functional which 
is extremized on the right-hand side of eq. (3.5) matches the standard expression for 
the horizon entropy of a black hole. While this proposal passes a variety of consistency 
tests, e.g., see [46, 47, 49], there is no general derivation of this holographic formula 
(3.5). However, a derivation was recently provided for the special case of a spherical 
entangling surface in [21]. 

In [23, 49], the above expression (3.5) was extended to holographic theories dual to 
GB gravity (2.9) in the bulk. The new prescription still extremizes over bulk surfaces 

V which connect to the entangling surface at the asymptotic boundary, however, the 
entropy functional to be extremized becomes^ 



S = ^ I d'^-'xVh 



1 H TZ 

(d-2)(rf-3) 

9, (d-2)(rf-3) ^ ^ 



Here, h [h) is the induced metric on (the boundary of) the bulk surface v, TZ is the 
Ricci scalar of this induced geometry and /C is the extrinsic curvature of the boundary 
dv at the asymptotic cut-off surface. Note that we only apply this expression for d > 4 
since it is only for these dimensions that the GB interaction (2.10) contributes to the 
gravitational equations of motion. Of course, if we set A = in the above expression. 



^We are using 'area' to denote the ((i-l)-diniensional volume of v. If eq. (3.5) is calculated in 
a Minkowski signature background, any extremal surfaces are saddle points of the area functional 
and one should choose the extremum with the minimum area. However, if one first Wick rotates to 
Euclidean signature, the extrcmization procedure actually corresponds to finding the global minimum 
of the area functional. 

^This expression was motivated by the construction of black hole entropy for Lovelock gravity 
appearing in [50]. We note that when evaluated on a general surface this functional will not match 
the Wald entropy [51]. However, the two agree when evaluated on a stationary black hole horizon. 
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it reduces to 27iA{v)/ip^^ and we recover eq. (3.5). Note that the extrinsic curvature 
term in eq. (3.6) plays the role of a 'Gibbons-Hawking' surface term to ensure that the 
variational principle is consistent. 




r = -txj y 



Figure 1: (Colour Online) The strip V in the asymptotic boundary, with the minimal surface 
V in the bulk ending on dV. The entangling surface dV consists of two flat (hyper)planes 
positioned at = i/2 and = —£/2. A regulator length H is introduced to limit the extent 
of these planes along the remaining directions. 

Now let us begin to consider evaluating the holographic entanglement entropy for a 
general RG flow. As in the previous section, we assume the bulk metric takes the form 
given in eq. (2.3). Then the boundary geometry is simply fiat space and we define the 
entanghng surfaces as follows: First recall that the entangling surface divides a Cauchy 
surface {e.g., the constant time slice, = t = 0) into two regions. As described above, 
we wish to consider an interval of length i and so we introduce two flat (hyper)planes 
at = i/2 and x^ = —i/2, as shown in figure 1. We also introduce a regulator length 
H to limit the size of the two planes along the x^'^' " ''^"^ directions, e.g., we can imagine 
the boundary is periodic in these directions with length H ^ i. Hence the area of either 
plane is H'^~^, as described at eq. (3.2). In calculating the holographic entanglement 
entropy, we consider bulk surfaces that end on the entangling surface as r — oo, as 
shown in figure 1. With the 'slab' geometry described here, the radial profile of these 
surfaces will only be a function of the coordinate x = x^. We will write the profile as 
r = r{x,w = i) where w indicates the width of the interval which sets the boundary 
condition, i.e., in the present case, x — )■ i/2 as r — )■ oo. Of course, the holographic 
calculations are only well-defined if we introduce an asymptotic cut-off surface as some 
r = rc. The position of this surface is related to a short distance cut-off in the boundary 
theory, i.e., Vc = L \og{L/6). The radial profile will define another useful UV scale e 
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with Tc = r{x = {i — e)/2,w = £), i.e., the profile intersects the cut-off surface r = at 
X = {£ — e)/2. Another useful scale in the bulk surface is the minimal radius which 
it reaches in the bulk, which appears as = r{x = 0,w = i). 

Given the background metric (2.3) and our ansatz r{x, t) for the profile of the bulk 
surface, we find that eq. (3.6) reduces to the following simple expression 



where r = dr/dx and A' = dA/dr. Now one may treat the above expression as an 
action which is varied to find a second-order differential equation to determine the 
profile r{x,i). However, since the integrand above has no explicit dependence on the 
coordinate x, the following is a conserved quantity along the radial profile^^ 



e 



-dA ^^2A _j_ p^3/2 



KM = e2^ + (l-2AL2/l'2)r2 ' ^^.S) 

This leaves us with a first-order equation for the profile, which should be easier to 
solve. In principle then, our goal is to solve for r{x,i) in a given holographic RG fiow 
geometry, i.e., for a specific conformal factor A{r), and then substitute the solution 
back into eq. (3.7) to calculate the entanglement entropy. 

Before going on to consider the entanglement entropy and c-function for RG fiow 
geometries, let us first examine the results when the bulk geometry is simply AdS 
space, i.e., at a fixed point of the fiow where the boundary theory is conformal. Recall 
that for the AdS vacuum A{r) = r/L. Let us begin by setting A = and considering 
the results for Einstein gravity in the bulk.^^ The case of three-dimensional AdS or 
a. d = 2 boundary CFT is special since the entanglement entropy yields a logarithmic 
UV divergence 

If we recall that the central charge of the boundary CFT is given by c = 127r L/£p, we see 
that this expression precisely reproduces the expected result (1.3) for the entanglement 
entropy of a two-dimensional CFT. Next turning to Einstein gravity with o? > 3, the 
entanglement entropy for the interval is given by [45] 

47r L'^-WhV-' 2''n('^+'y' fn^))Y~' L''~WhV-' 



~ d-2 \ b) id-2) l^r(^^ I ^'-^ • ^^-^^^ 

^''if we denote the integrand in eq. (3.7) as C, then Kd{i)^^ = C — ^ f. 
^^Note that in this case, the AdS curvature is given by simply L = L. 
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Here we see the general structure given in eq. (3.2) with two terms, a power law diver- 
gence proportional to {H/5Y~'^ and a finite contribution proportional to {H/lY~'^. 
Next for A = 0, both of the central charges in eqs. (2.13) and (2.14) are identi- 
cal and we use this fact to define Cd in eq. (3.2) for Einstein gravity: Cd{\ = 0) = 



7r'^/Vr(rf/2) (L/4 



\d-l 



T 



As described previously then, we can extract this 



central charge from the above entanglement entropy using eq. (3.3), which yields 



l)d-l 



Cd = l^d 



OS, 



1 



d-1 



CFT 



i7"'-2 di 



with 13d 



^2(d-l)' 



7r2'ir(d/2) \ T( 



d 



(3.11) 



2(d-l)' 



Hence we have identified the precise value of f3d (for d > 3) which appears as the 
coefficient in eq. (3.4) of the c-function. 

Finally let us apply the above formulae to calculate holographic entanglement en- 
tropy with the strip geometry for the boundary CFT dual to the AdS vacuum in GB 
gravity (2.9). To simplify the final results, it is convenient to first treat r as the 
independent variable, in which case to fix the profile of the bulk surface, we must de- 
termine x(r). Next we choose a new radial coordinate r = e"^/^ and define 
h = e^^/^ (drxY- Note that as r — )• oo, r — t- and further one can show at r = rm, 
r = 1. Now with these choices, eq. (3.8) becomes 



r 



d-1 



1 + hf/^ 



v^(l-2A/oo + /i) 



1. 



(3.12) 



In general, this equation yields three roots for /i(r) and the relevant solution is the 
real root which can be continuously connected to the A = solution: h = T'^^^~'^y{l — 
^2{d-2)y Now it is straightforward to see that the entanglement entropy (3.7) can be 
written as 



5, 



d-l 



d-2 



CFT 



d 



(1 + 2A/, 



X I 



d-2 




d-2 



l + {d-2) 



(3.13) 



where / = dry/h. Then applying eq. (3.11), we can express the central charge in 
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the finite contribution as 



12 



Cd{X) = 47r/34l + 2A/, 



rd-l 

(jd-l 
--P 



rd-2 



1 + (ci - 2) 



dr 



1 + 2A/oo + h 
(1 + 2A/oo) VT+ 




Regrettably, we do not have a closed analytic expression for Cd{\) in terms of the two 
central charges Ct and a^. Hence we have numerically evaluated the above expression 
and plotted Cd{X)/a^ as a function of CT/a*d in figure 2a for several values of d. Note 
that in this figure, Crf(A)/a^ = 1 at Ct/cl"^ = 1 for all of the values of d since this 
corresponds to A = or Einstein gravity in the bulk. From these curves, we can 
infer that Cd{X) is a complicated nonlinear function of both Ct and ajj. We can also 
illustrate this fact as follows: In the vicinity of A ~ or Ct 
linearized analysis of eq. (3.14) to find 



a^, we can make a 



Crf(A) = Crf(O) 1 



3{d-l] 



A + 0(A^ 



where 



Cl(A) 



{d-l)id-2) 



C 



T 



d{d - 3) 



(3.15) 



(3.16) 



2 ^ 2 

Here, we have defined Cl(A) as the linear combination of the two central charges in 
eqs. (2.13) and (2.14) which yields an 0(A) expansion which precisely matches that for 
Crf(A). Next, we consider the ratio of Crf(A) and Cl{X) over the full (physical) range of 
A. Since Cl{X) can vanish in this range, it is convenient plot the ratio Cl{X) /Cd{X) as a 
function of CtI^^ as shown in figure 2b. This figure illustrates even more dramatically 
our previous observation that Cd(A) is a complicated nonlinear function of both Ct and 
a^. At this point, let us add that since Cd{X) 7^ ajj, the central charge identified in [20] 
as satisfying a c-theorem, we might not expect that our new effective central charge 
C(i(A) will always flow monotonically in holographic RG flows for general A. 



4. Holographic flow of c-function with Einstein gravity 

In this section, we examine the behavior of the c-function (3.4) in a general holographic 
RG flow dual to Einstein gravity. We will first discuss the flow of the c-function in d = 2 
and then generalize it to arbitrary dimensions. 

^^Note that analogous results were given for the case d = 4 in [23]. However, we note that the 
calculations presented there did not include the 'Gibbons-Hawking' surface term in eq. (3.6) and 
hence their expressions do not match those presented here. However, we have verified numerically 
that the effective central charge in [23] agrees with eq. (3.f4) when d = 4. We also observe that the 
leading divergent term in eq. (3.f3) is proportional to Ct while without the 'Gibbons- Hawking' term, 
this term is proportional to a*^. 
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Figure 2: (Colour Online) Panel (a) plots Cd{X)/a*^ as a function of Ct/cl*^ for GB gravity, 
while panel (b) is a plot of Cl{X) /CdiX) as a function of Ct/o,*^- Note that both Cii{\)/a*^ = 
1 = Cl{X) /Cci{X) at Ct/cl'^ = 1 which corresponds to Einstein gravity in the bulk {i.e., A = 0). 
Each curve runs over the physically allowed range of Ct/cl^ for the given value of d — see 
eq. (2.15) for more details. 



For d = 2, the bulk theory is three-dimensional Einstein gravity coupled to, e.g., 
a scalar field with a nontrivial potential, as described in section 2. Our holographic 
expression (3.7) for the entanglement entropy of a strip can be written as 



S 



4:71 

JO 



dx \l\ 



j.2 _|_ g: 



2A 



(4.1) 



Note that the integration above runs over half of the range, i.e., x G [0, (i* — e)/2]. 
Further the conserved charge (3.8) simplifies to 



,2A 



(4.2) 



To calculate dS/di, we note that r = r(x, i) and A = A{r{x, i)), i.e., the profile of the 
extremal surface implicitly depends on the strip width i. If we vary S with respect to 
£, keeping the UV cut-off r = Vc fixed, we will get two contributions: one coming from 
change in the limits of the integration and second from change in the solution r{x,i). 
We write them as 



dS 



An 



1 



di ip ^^2 + g2A(r) 

where we have used the equation of motion for r{x,i) 

r-A' (2r2 + ^2A^^ ^ g 



dr 

m 



(4.3) 



(4.4) 
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to cancel the bulk contribution. Since the UV cut-off is fixed while performing the 
variation, we get some extra constraints between f = dxV and r' = dir at the asymptotic 
boundary. Taking variation of relation r(^,£) = Tc with respect to we get 



rue. 



1 - 



de 
di 



dr{x, tj 
di 



0. 



(4.5) 



Substituting this relation, as well as eq. (4.2), into eq. (4.3) gives us following expression 
for dS/di: 

dS _ _ 4:71 Idr 
~dl ~ ~ipK2{i)^di 



(4.6) 



In the above relation, the partial derivatives of r(x, i) are evaluated near the asymptotic 
boundary. To further simplify this expression, we use the Fefferman-Graham expansion 
near the boundary [30]. In terms of the radial coordinate r, this expansion takes the 
form [24] 



ds' 



dr' 



^^'^^ f{r)r],jdx' dx^ 



(4.7) 



where 



/(r) = 1 + ai0^ e-'""/^ + 03^^ e"^"^^/^ + ■ ■ ■ . (4.8) 

In this expansion, L is the AdS radius in the UV region {i.e., as r — )■ 00) and a = rf — A 
where A is the conformal weight of the operator dual to the bulk scalar field. Near the 
boundary, the coordinate r is very large and hence it is sufficient to work with only 
the leading order term in the expansion (4.8). Although A{r) has a complicated profile 
deep inside the bulk, near the boundary it will have the simple form 

A{r) = r/L. (4.9) 

For this A{r), eq. (4.2) can be re-expressed as the following equation of motion: dx/dr = 



e '^^/^ / \ K2 — e~'^^'/^. The latter is easily integrated to yield the following solution 



X 



e-2r/L _ 



(4.10) 



where the integration constant was chosen so that x— )-£/2asr— j-oo. Next we 
differentiate the above solution with respect to x and I to find f and dr/dl, treating 
that K2{i) as a function of i - see appendix A for further details. Taking the limit 
r — )■ 00 in ratio of dr/d£ and f appearing in eq. (4.6), we find that 



1 dr 
r di 



x=£/2 



(4.11) 
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This relation not only simplifies eq. (4.6) but also ensures that the first derivative of S 
is indeed finite for all RG flow solutions. Using this relation in eq. (4.6), we arrive at 
following elegant form of the c-function (3.1) for arbitrary RG flow backgrounds: 



C2 = 3i 



dS 



67r £ 



(4.12) 



The next step is to show that this c-function increases monotonically along holo- 
graphic RG flows. Implicitly, the extremal bulk surfaces v on which we are evaluating 
the Ryu-Takayanagi formula (3.5) extend to infinite r at a; = ib£/2 and pass through a 
minimum r = at x = 0. The latter radius gives us an indication of which degrees of 
freedom the entanglement entropy is probing, i.e., for smaller values of r^, we expect 
the entropy and C2 responds more to the IR structure of the RG flow. Hence in the 
following, we will study behavior of C2 as a function of the turning point radius 
and we wish to establish the 'c-theorem' as dc^jdrm > - at least for background 
geometries that satisfy appropriate constraints. 

Comparing to the field theory construction of [2] , we note that there the c-theorem 
was formulated as dc^jdl < 0. Naively, this result matches with the holographic 
inequality which we wish to establish since we expect that as the width of the strip 
increases, the minimal area surface will explore deeper regions in the bulk geometry. 
The two inequalities would be rigorously connected if we could prove a second inequality 
di/drm < for consistent holographic models. However, as we will see in the next 
section, in fact this inequality does not hold for all extremal surfaces. However, we will 
still find dc2/di < in all cases of interest. The violations of the previous inequality 
are associated with unstable saddle-points which do not contribute to the physical 
entanglement entropy. Hence, in section 5, we will find that the behaviour of the 
entanglement entropy in general holographic RG flows provides a richer story than 
might have been naively anticipated. 

Returning to the flow of the c-function, we note that at the minimum of the bulk 
surface, we will have r{0,i) = and r{0,£) = 0. Hence considering eq. (4.2) at this 
turning point, we find 



Here it is natural to treat this constant of the motion as a function of r^, rather than 
i. We will also work with width of the strip i as function of r^- Then combining 
eqs. (4.12) and (4.13) yields 



= e 



(4.13) 




(4.14) 
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Now to express ^ in terms of r^, we begin with the relation 



dx = —= dr , =—. 4.15 



r 



Here in the final expression we have used eq. (4.2). Now above, we will apply integration 
by parts using 



2 _ e-2A 



2 



dr -^=^^ = JkI - . (4.16) 



to find that 

A" 



oo 



i = 2LK2 + 2 / dr—^Kl - e-^^ , (4.17) 

where we have used eq. (4.9) to evaluate A'{r) at r = oo. Further we can differentiate 
this expression with respect to to get 

^ff /-oo A" 1^2 

-2LA'{rr^)K,-2A'{r^) I dr — ^=^==. (4.18) 



dr^ -^-'l^^'A-^ 
Now substituting eqs. (4.17) and (4.18) into eq. (4.14), we find 



e~ 



dc^ 12 IT A' (rm) r , A" 

dr- 



drm 4i^2(r™) A'^^/Kf 



'2 ^-2A 



e 



^'^^'<-' ^'</.4;. (4.18) 



ipK2{rm) In A 

12 iTA'(r^) r , 1 



K2{rm) 



I dxjr2 {T\-T\)>0. 



In the second line, we have used eq. (4.2) to convert the integration over r to one over 
X. In the last line, we have used Einstein's equations to replace A" by the components 
of the stress tensor. As for the discussion of holographic c-theorems in section 2, the 
final inequality assumes that the bulk matter fields driving the holographic RG fiow 
satisfy the null energy condition. The latter ensures that the integrand is negative. 
The overall inequality also requires K2{rm) > and A'{rm) > 0. The first condition 
is obvious from eq. (4.13) while the second can be established as follows: Given the 
null energy condition, it follows that A" < which means that A' is everywhere a 
decreasing function of radial coordinate r. Implicitly, we are assuming the bulk ge- 
ometry approaches AdS space asymptotically, i.e., the dual field theory approaches a 
conformal fixed point in the UV. Hence with r — j- oo, we see the minimal value of A' is 
A' = 1/L, where L is the asymptotic AdS scale. Since this minimal value is positive, it 
must be that A' is positive everywhere along the holographic RG fiow. Hence dc2/drm 
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is positive and our two-dimensional c-function increases monotonically along the RG 
flow if the bulk matter satisfies the null energy condition. 

We now turn to proving the monotonic fiow of the c-function (3.4) for higher di- 
mensions. The required analysis is a straightforward extension of the above calculations 
with d = 2. In particular, one finds that eq. (4.12) generalizes to 



Cd 



2n/3d t 



d-i 



with d boundary dimensions. The conserved quantity (3.8) is now given by 



r2 



□ 2A(r) 



(4.20) 



(4.21) 



We have relegated the detailed derivation of eq. (4.20) to appendix A. However, we can 
see from this result that all the complexities of determining the c-function boil down 
to evaluating the conserved charge (4.21) for the minimal area surface. We might note 
that we can evaluate this expression at the minimal radius (where r = 0) to find 



KJr. 



dV'm 



-(d-l)A(r™) 



(4.22) 



which generalizes eq. (4.13) to general d. 

Combining eqs. (4.20) and (4.22), we further find 



dcd ^ 2id-l)7i (3di' 

drm Kr, 



d-2 



di 

d/Trn 



(4.23) 



To express i in terms of r^, eq. (4.15) now becomes 



dx 



dr 
r 



dr 



_ e-2{d~i)A A' 



(4.24) 



where the last expression follows by using eq. (4.21). To make further progress, we 
observe that 



dr 



A' e-'^^('-) 



-dA 



-2{d-l)A 



dK, 



d 



3d 



^-i,d-2)A j^2 



-2{d-l)A 



2' 2(d- 1)' 2(rf- 1) 

-{d-2,) A j^2 



-2{d-l)A 



-2{d-l)Air) 

d-2 



d 



(4.25) 

g_2(d-l)A(r)- 



2' 2{d-l)' 2{d-l)' 



We have presented the second expression above to illustrate that this result is simply an 
extension of eq. (4.16) for general d but in the following, we will use the more compact 
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expression given in the first line. With eq. (4.25), we can integrate by parts in eq. (4.24) 
to find 



■id-2 
2(d-l) 



dA'hr 



r 



2d~l 



(4.26) 



2(d-l)^ 

°° A" e~'^^ 



d 



3d 



2' 2(rf- 1)' 2(rf- 1)' 



Differentiating this result with respect to and making various simplifications yields 
di K'/'-'^ ^ {w^) 



dVr. 



d 



+2A'{rm) 



2d-l 



2A'(r„ 



dr 



2{d~l) ^ 

A" e-'^^ 



-dA 



dv - 

„, A'^ JIC - e-^C'^-i)^ 



(4.27) 



d 



3d 



dKd [2' 2(rf- 1)' 2(rf- 1) 
Now substituting eqs. (4.26) and (4.27) into eq. (4.23), we find 



-2{d-l)A{r) 



dcd 
dr„ 



An{d-l)Pdr-'A'{rm) 



p -^di'^m) 



OO j^// 

dr 



^dA 



A7r{d-l)l3dr-'A'{rm) 



•d-2 A/ 



dx 



A^ 



2 _ g-2((i-l)A 



(4.28) 



4 7r/3dr-M'(r 



dx — {T\ - T\.) > . 



The steps here are essentially the same as in our analysis of eq. (4.19) with d = 2. The 
key requirement for the final inequality to hold is that the bulk matter fields driving the 
holographic RG fiow must satisfy the null energy condition. With this assumption then, 
dcd/drm is positive and our (i-dimensional c-f unction increases monotonically along the 
RG fiow for holographic boundary theories dual to Einstein gravity in the bulk. 



5. Explicit geometries and Phase transitions 

In this section, we consider some simple bulk geometries describing holographic RG 
fiows. This allows us to explicitly demonstrate that the c-function (3.4) indeed fiows 
monotonically for boundary field theories dual to Einstein gravity. However, we will also 
find that for some RG fiows, there is a 'first order phase transition' in the entanglement 
entropy as the width of the strip i passes through a critical value. Technically, denoting 
the behaviour in the entanglement entropy as a phase transition is inappropriate - after 
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all, the system itself, i.e., the state of the boundary field theory, does not change at 
all. However, as we will see below, in our holographic calculation of the entanglement 
entropy, there are competing saddle points and the dominant saddle point shifts at a 
critical value of the width. Of course, this behaviour is reminiscent of that seen in 
holographic calculations describing thermodynamic phase transitions [52] and so we 
adopt the nomenclature 'phase transition' to convey this picture. The phase transition 
is first order and so the entanglement entropy is continuous at the critical width if, 
however, the derivative dS/di is discontinuous at this point. As a result, the c-function 
drops discontinuously at the phase transition. 

Implicitly, we are assuming that the holographic RG flows studied below are solu- 
tions of Einstein gravity and hence the entanglement entropy is determined by eq. (3.5). 
Explicitly, our RG flow geometries take the form given in eq. (2.3) and so are defined 
by giving the conformal factor A{r). Here we note that in all of the examples we 
consider, A"{r) < and so the geometry could solve Einstein's equations with matter 
fields satisfying the null energy condition. In appendix C, we consider one approach 
to constructing an appropriate scalar field theory that could realize the latter. In any 
event with A"{r) < 0, the holographic c-theorem of section 4 will be satisfied. That 
is, dcd/drm > or alternatively, the c-function decreases monotonically as the corre- 
sponding extremal surface extends deeper into the bulk geometry, as will be shown 
below. 

In general, we will consider arbitrary values of the boundary dimension in the 
following. However, to begin, we consider a very simple example of a step flow and 
the discussion will be limited to the case d = 2, i.e., a three-dimensional bulk. The 
step profile consists of two AdS geometries with different curvature scales are patched 
together at some finite radius. With such a simple profile, the behaviour of the entan- 
glement entropy and the c-function can be determined analytically. Our analysis with 
c? = 2 is easily extended to higher d but we do not present the results here. In the 
subsequent subsection, we also examine smooth profiles describing an holographic RG 
flow and allow for arbitrary d. However, numerical analysis is required to understand 
the behaviour of the entanglement entropy for these smooth profiles. 

5.1 Step profile 

We limit the discussion here to three-dimensional gravity and consider a bulk geome- 
try^^ which patches together two AdS regions with different curvatures at some finite 
radius r = tq. Using the metric ansatz in eq. (2.3), the conformal factor A{r) is given 

^■^Various aspects of the flow of entanglement entropy in this example was also studied in [25] for 
= 2, 3 and 4. 
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by 

. . ^ r AiMr) = ^ + for r < ro 

where Lj/y and L/^ correspond to the AdS radius in the UV and IR regions, which 
we denote as AdSjyy and AdS/R in the following. The constant term added to A{r) 
in the IR region ensures that the conformal factor is continuous at r = tq. Of course, 
it is not differentiable there and some stress energy with 5-function support would be 
required to make this geometry a solution of Einstein's equations. As discussed in 
previous sections, there is a conserved quantity (4.2) which plays an important role 
in determining to the entanglement entropy and the c-function. Clearly, there are two 
classes of minimal area surfaces in this geometry, namely those that stay only in AdSiyy 
and those that penetrate deep enough into the bulk so that the minimal radius is 
in AdS/ij. In either case, the conserved quantity K2 is given by eq. (4.13) and hence 
we have Kuy = e~"^^^^^"^^ for > tq and Kjji = e~'^'^^'^^'^ for < tq. To regulate 
the entanglement entropy, all of these surfaces are terminated at a large cut-off radius 
r = Tc in the UV region. 

To find the minimal area surface, we will solve eq. (4.2) for x = x{r). First, we can 
invert the latter equation to find 

^ = , =. (5.2) 

dr ^Kl - e-2^(^) 

Above we have discarded the root with an overall minus sign because we will only 
consider the branch of solutions covering the interval x e [0, £/2] in the following, for 
which dx/dr > 0. For the minimal surfaces that stay entirely in AdS[/v, we can easily 
integrate (5.2) to find 



X = Luv^Ky-e-^-/^uv, (5.3) 

using Auv as given in eq. (5.1). The integration constant is chosen here so that x = 
at r = Tm, which also implies that 

i = 2LuvKuv = 2Lc7ye-'''"/^^^ . (5.4) 

As noted above, this solution is valid for > vq, which implies i < £2 where 

£2 = 2Lc7ve-'^»/^^^ . (5.5) 

Next we turn to the second class of minimal area surfaces, which penetrate into 
AdS/R. In this case, we have K2 = Kjr = e'"^'^^"^"^^ with the turning point in AdS/ij. 
We divide the relevant solutions of eq. (5.2) in two parts: xmir) G [0,Xt] describes the 
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portion of the extremal surface in AdS//j and Xuvi^^) ^ l^t,^/'^] represents the part in 
AdSf/y. Here we have defined the transition point Xt such that Xt = xjR^ro) = xuvi^o). 
Now integrating eq. (5.2) with the appropriate conformal factor (5.1) for each segment, 
we find 



XiR = L,HA/if|^-e-2A.«M 

i 
2 



Xuv = LuvJk^ir - e-2^^vM + _ _ LuvKir . (5.6) 



Above, the integration constants were chosen so that xjr = at r = and xuv = 1/2 
as r — )■ oo. Combining these solutions at x/ij,(ro) = xuviro) = xt, we find 

Xt = J ^^^j — [LuvKiR - ^ ) , (5.7) 



Luv — LiR \ 2 ^ 

e = 2LuvKjR - 2{Luv - Ljr) ^ K]^ - e-^^^/^uv . (5.8) 

Here we have ensured that the solution (5.6) is continuous at r = ro but the first 
derivative is also continuous at this transition point because of the form of eq. (5.2) 
and the continuity of the conformal factor. Implicitly, eq. (5.8) gives the relation 
between i and since Kjr = e~^^«('''"). As the physically relevant quantity in the 
boundary theory is i, we invert this relation to find the following two solutions for Kjr: 

Luvi ± (Luv - Ljr)^P - ALjr{2Luv - LjR)e-^^o/Luv 

KiR± = , . (5.9) 

'^^^my^-i^uv — ^ir) 

Above, both of these roots provide real solutions for i > ii with 

£1 = 2^Lir{2Luv-Ljr) e-^o/^^^ . (5.10) 



It is also useful to define ri, the value of the minimum radius at £ = £1, for which we 
find 

V Lir{2Luv — Ljr) 

Now the root Kjr^ is a monotonically increasing function of i over the range ii < i < 
00 and for any i in this range, there is a consistent solution for the extremal surface. 
The corresponding values of the minimal radius are ri > > —00. Now the second 
root KiR- decreases for ^ ~ £1, however, it is an increasing function for large values of 
i. Kjr- has a single minimum at i = £2, i.e., precisely the width defined in eq. (5.5) 
for the discussion of solutions remaining entirely in AdSj/y. At this minimum, Kjr- 
takes the value e~^'^^'^°^ = e"^"/^*^^. We find that Kjr_ yields a consistent solution 
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for the extremal surface as long as £i < i < £2- However, for £ > £2, the solutions 
corresponding to Kjr- are inconsistent, e.g., dx/dr is not positive throughout the 
range x G [0, i/2]. We note that for the consistent solutions, while £ runs from £1 to £2, 
the minimum radius of these surfaces increases from ri to tq. That is, in contrast 
to the previous two families of solutions, here we have d£/drm > 0! 

Hence the following picture has emerged for the extremal surfaces: Beginning with 
the minimal radius in the range 00 > Vm > tq, there is a family of extremal solutions 
which remain entirely in AdSjyy. As can be seen from eq. (5.4), the width £ increases 
monotonically as decreases, reaching the maximum £ = £2 when = tq. Below this 
point, we make a transition to a new family of solutions which begin to penetrate into 
AdS/ij. For tq > rm > Ti, the relevant family of extremal surfaces corresponds to the 
branch with Km-- In this regime, £ actually decreases as continues to decrease, i.e., 
d£/drm > 0. When reaches ri, as given in eq. (5.11), £ = £1 and we make another 
transition to the third family of extremal surfaces. These solutions correspond to the 
branch with i^/_R+. In this regime ri > > —oo,£ again increases monotonically 
as Tm decreases. Figure 3 illustrates this behaviour for all three families of extremal 
surfaces. Now for any particular value of the turning point radius r^, we see there is 
unique extremal surface. However, if we consider the solutions as a function of the strip 
width £, there is a unique solution for £ < £1 and £ > £2. In the intermediate range 
£1 < £ < £2, there are in fact three possible extremal surfaces for any given width. 
Given three possible saddle points, we are instructed in eq. (3.5) to find the extremal 
surface with the minimum area in order to evaluate the entanglement entropy. This 
situation with multiple saddle points is also the typical scenario that one encounters 
in the holographic description of a thermodynamic phase transition [52] and in fact, 
we will find the latter extends to the present situation. That is, we see below that the 
entanglement entropy undergoes a 'first order phase transition'. 

Hence having found the solutions for the extremal surfaces, we will present an 
entropy, i.e., S = 27rA{v)/£p, for each of these surfaces. But, of course, in the interme- 
diate regime described above, the true entanglement entropy is given by the solution 
which minimizes this quantity. Let us begin with the solutions (5.3) which remain 
entirely in AdSjjy. For this case, the entropy turns out to be 



Here the result is expressed in terms of a short-distance cut-off in the boundary theory 




(5.12) 
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Figure 3: (Colour Online) This plot illustrates the behaviour of £ as a function of for all 
three families of extremal surfaces. The red, blue and green portions of the curve correspond to 
the contributions coming from S'//j(ii'/R+), SjR{Kjji^) and Suvi respectively. In particular, 
we see that dl/drm < for Sir{Kirj^) and Sjjv-, while M/drm > for 5/r(-R'/r_). For this 
plot, we chose Ljjv = 1; Ljr = 0.1 and ro = 0. 



6, which is related to the radial cut-off by^'' 

Tc = Luv log (Luv/S) . 



(5.13) 



Now let us consider the extremal surfaces given by eqs. (5.6) to (5.9), which penetrate 
into AdS/ij. For this case, the calculation of entanglement entropy results 
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(5.14) 



-2xt 



We convert the cut-off e above to 5 using r(x = (£ — e)/2, £) = r^, which yields 

5^ 



LuvKiR 



(5.15) 



^""^This matches the standard cut-off z„„;„ = 5 in Poincare coordinates where the AdS metric takes 
the form ds"^ = {Lfjy/z'^) {r]i.jdx^dx^ + dz^) ~ e.g., see [21]. 
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Using this relation, we can write 
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(5.16) 



This result is valid for both of the roots Kir± given in eq. (5.9). 

Now as described above, the UV family of solutions (5.3) provide the unique ex- 
tremal surface for any £ < ii and hence the entanglement entropy is given by 5* = Suv 
in this regime. Similarly, for £ > £2, the extremal surface is again unique and hence 
the entanglement entropy is given by S" = Sir{Kir^). In the intermediate regime 
£1 < £ < £2, have three extremal surfaces and we must identify which of these yields 
the minimal entropy. In particular, we always find SiR{Kfji+) < Sir{Kir^) and hence 
the branch with Kjji_ never plays a role in determining the physical entanglement en- 
tropy. Therefore the latter is found by comparing Suv and Sir{Kir+). It turns out 
that these two entropies are equal for some critical width £t with £1 < £t <i £2- Further 
Suv > Sir{Kjr+) for £ > it and Suv < Sir{Kiji^) for £ < £t. Hence we find that, the 
entanglement entropy for the step profile (5.1) is given by 



S 



Suvii) for £<£t 

SiR{Kin+m ioT£>£t, 



(5.17) 



with £t defined by Suvi^-t) = SiR{Kjfi^{£t)). In particular, we observe that the en- 
tanglement entropy exhibits a first order phase transition at the critical width £ = £t. 
Further, we note that at transition, S{£) is continuous but not differentiable. We have 
illustrated all of this behaviour in figure 4, which plots S{£) — S{£2) versus log(£/£2) for 
specific values of the parameters, Luv, Lj^ and rg, defining the profile. 

Given the entanglement entropy (5.17), we turn to the calculation of the c-function 
defined in eq. (3.1). We could proceed here by explicitly differentiating the various 
expressions above, e.g., eq. (5.16), with respect to £ to determine C2. However, this 
calculation only verifies the final result which was already determined in our general 
analysis in section 4, namely eq. (4.12). In fact, this result applies for all three families 
of extremal surfaces and so we have 



C2 



&K_e_ 
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(5.18) 



Kin^ from Sir{Kir+) , 
"'^^Since S{£) diverges as 5 — > 0, wc plot the difference S{£) — S{£2) which is independent of S. 
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Figure 4: (Colour Online) This plot illustrates the entropy for all three families of extremal 
surfaces as a function of the width of the strip i. The green curve corresponds to Suv', the 

red, to Sifi{Kffi and the blue, to Sm^Kfji^). The phase transition occurs at i = it 

where the red and green curves cross. For this plot, we have chosen Ljjv = !> -^/fl = 0.5 and 
ro = 0. 

where Kjjv = and Kjji± are given in eq. (5.9). Figure 5a plots C2 as a 

function of the turning point radius - or rather e^™/^'^^. We see that dc2/drm > 
everywhere in the figure, which is again in keeping with the expectations of our general 
analysis in section 4. However, because of the phase transition, not all values of are 
relevant for the c-function (3.1) evaluated on the physical entanglement entropy (5.17). 
In figure 5a, the region between the vertical dashed lines is excluded and the physical c- 
function jumps discontinuously between the values at the points labeled A and B. This 
behaviour is also illustrated in figure 5b where the c-function in eq. (5.18) is plotted as 
a function of the ratio i/if The phase transition ai i = it again takes C2 between the 
points labeled A and B, which now lie on the same vertical dashed line in this figure. 
That is, at this critical value of the strip width, the c-function drops from the value 
given by Suv {on the green curve) to that given by Sifi{Kjfi^) (on the red curve). 
Again this discontinuity arises because of the first order nature of the phase transition, 
i.e., the entanglement entropy is continuous but not different iable at this point. If we 
consider only the physical values of C2, then we also find dc2/di < in keeping with the 
general expectations of field theory analysis of [2]. Of course, figure 5b also illustrates 
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Figure 5: (Colour Online) Panel (a) plots the c-function (5.18) as a function of e~'^^l^uv _ 
This plot clearly illustrates that ci decreases monotonically as decreases, in accord with 
the analysis of section 4. However, the region between the vertical dashed lines is excluded 
by the phase transition in the entanglement entropy. Rather the physical c-function jumps 
from the point labeled A to that labeled B. Panel (b) plots the c-function with respect to 
ijlt- One can see that C2 is multi-valued in the region i\ < i < £2, as noted in eq. (5.18). 
Again &t I = li^ the c-function drops discontinuously between the points labeled A and B. 
The red, blue and green portions of both curves correspond to the contributions coming from 
Sjr{Kjji^), SjFi{Kjji^) and Suv, respectively. For this plot, we chose Luv = Ij Ljr = 0.1 
and ro = 0. 



that dc2/di > on the branch associated with Km-. However, as emphasized above, 
this family of saddle points is not relevant of the physical entanglement entropy (5.17). 
The 'unusual' behaviour of the c-function on this branch arises because M/dr^-^ > 
for this family of solutions. Given the behaviour illustrated by this simple example, it 
seems likely that in general any branch of extremal surfaces with the latter property 
will correspond to unstable saddle points which are not physically relevant. 

5.2 Smooth profiles 

The simple example in the previous section has alerted us to the possibility that the 
entanglement entropy S{i) may experience a phase transition with respect to changing 
the strip width i. However, one should worry that this result is an artifact of the 
artificial shape of the step profile in eq. (5.1). Hence we consider some smooth profiles 
in this section and examine to what extent this phase transition survives for these more 
realistic holographic RG flows. Again our definition of the holographic entanglement 
entropy is given in eq. (3.5) and so implicitly we are assuming that the bulk geometry 
is a solution of Einstein's equations. In appendix C, we consider the scalar field theory 
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that would be necessary to realize the latter. In this section, we will consider arbitrary 
values of the boundary dimension d. 

Let us first consider a smooth flow between the UV and IR fixed points with the 
following conformal factor 



(2 cosh(r/i?) ) ^, (5.19) 



Notice that A{r) ~ r/L — ^rjR in the limit r — )■ +cxd and Air) ~ rjL + ^rjR for 
r — i- — oo. Hence the geometry approaches AdS space in both of these limits with 

1 _ 1 7 1 _ 1 ^ 7 ^^2^^ 

L R Ljji L R 

The parameter R controls the sharpness of the transition in the holographic flow be- 
tween the UV and IR fixed points while the change in the AdS scale is controlled by 
the combination j/R. In the limit i? — with 7/-R fixed, we would recover a step 
profile of the form given in eq. (5.1). 

To proceed further in examining the possibility of a phase transition, we used the 
analysis presented in previous sections and examined the extremal surfaces numerically 
for the above holographic flow profile. First, using eq. (4.21), the equation determining 
the shape of the extremal surfaces is reduced to a first order equation, 

^ = , (5.21) 

as appears in eq. (5.2) for d = 2. Then families of surfaces are easily constructed as a 
function of the turning point radius using eq. (4.22). Numerically integrating from 
the turning point out to the asymptotic region, we can then determine i{r„i). 

Let us add a few more details about the numerical analysis: Near r = r^, one finds 
that X ~ y/r — and hence eq. (5.21) is singular precisely at r = r^, the putative 
starting point of our numerical integration. So to simplify the numerical analysis, we 
define 

y{r) = y/r - rmx{r) , (5.22) 
for which the equation of motion (5.21) becomes 

dy ^ y y-r - e"'^^^^^ 

dr 2(r-r„,) ^Rj - 6-2(^-1)^(0 ' ^ ' ^ 

With this new coordinate, y ^ {r — Tm) near r = rm and the right hand side of the 
equation of motion (5.23) is finite. To set the initial conditions, we use eq. (5.23) to 
find the leading terms in a series expansion of y{r) in r — r^- Now we can numerically 
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Figure 6: (Colour Online) Plot of strip width £ as a function of e^^l^vv fQj- various values 
of the profile width K. With small i?, there is a regime where dljdrm > 0. However, for large 
R, M/drm < everywhere. There is a critical value Rc for which di/drm reaches zero at a 
single point and is otherwise negative. The plot was prepared using the smooth profile (5.19) 
with L = 0.66 and 7 = 0.5i?, as well as d = 2. These parameter values yield Ljji = 0.5 and 
Luv = 1- 

integrate eq. (5.23) out from the turning point r = to large asymptotic values of r 
and find the the strip width £ using the relation 



Now as discussed previously, the appearance of a phase transition is directly re- 
lated to the appearance of a regime where dl/dr^ > 0. Figure 6 illustrates that such 
behaviour still arises for a range of parameters in the smooth profile (5.19). However, 
as shown in the figure when R grows (holding 7/-R fixed), this region decreases and 
eventually di/drm < for all values of r^. That is, there exists a critical value Rc such 
that, for R < Rc there is a first order phase transition while for R > Rc, we observe 
a smooth cross-over. At precisely R = Rc, there is a single point where di/dr^ = 
and the slope is otherwise negative. In this case, the phase transition would be second 
order. 



lim 




(5.24) 



m 
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Of course, one can go beyond the above analysis to identify the precise point where 
the phase transition occurs given specific values of the parameters in eq. (5.19) which 
produce a regime where M/drm > 0. Hence for some range of the width i, there will 
be multiple surfaces which locally extremize the entropy functional (A.l). Determin- 
ing which surface provides the dominant saddle point requires carefully regulating the 
entanglement entropy and comparing the values of finite parts of entropy for the com- 
peting saddle points. This analysis is essentially the same as in section 5.1, however, 
in the present case with a smooth conformal factor, the profile x{r) and the entropy 
integral are evaluated numerically. We will not present any of these results here. 

Let us now turn to the asymptotic expansion of the profile (5.19) and compare it 
to the Fefferman- Graham (FG) expansion given in eqs. (4.7) and (4.8). We find that 

/(r) = (1 + e-'^/^) ^ 1 - 7e-2^/^ + . . . . (5.25) 

which yields a = Luv / R in eq. (4.8). Now as described in section 2, the natural 
holographic interpretation of this flow would be that the UV fixed point is perturbed 
by a relevant operator, which would be dual to by a scalar field in the bulk theory. 
However, as discussed in appendix C, applying this interpretation to the bulk solution 
yields an upper bound on the parameter a appearing in the FG expansion (4.8), i.e., 
a < I + 1. This bound then becomes a constraint on i?, the width of the holographic 
profile (5.19). That is, 

(5,2«) 

Previously we found that the smooth holographic RG flows described by eq. (5.19) 
will still yield a first order phase transition in the entanglement entropy provided the 
width is sufficiently small, i.e., R < Rc- Hence the lower bound given in eq. (5.26) 
creates a certain tension. Namely, if Rc does not satisfy this lower bound, it seems 
likely that the phase transition is still an artifact of the artificial shape of the profile in 
eq. (5.19). We examine this question in figure 7a, where we have plotted {d+2)Rc/2Luv 
for different values of Luy /Li^. The bound (5.26) implies that [d + 2)R/2Ljjy > 1 
and as the figure illustrates, the latter can only be satisfied for sufficiently large d, i.e., 
d > 6. Hence the possibility of a phase transition is called into question for the physical 
dimensions d = 2,3, 4. 

Now the profile (5.19) was constructed to give a simple example which would 
smooth out the step potential studied in the previous subsection. One can easily 
generalize this profile to include more independent parameters and study the effect 
on the phase transition. Hence as another simple example, we consider the following 
conformal factor 

^Air) ^ ^r/L ^^2r/R ^ 3^ + e'^^/^) . (5.27) 
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Figure 7: (Colour Online) Panel (a) plots ^ 11^''^'' ^ function of 4^ for various dimensions 



{d+2)R. 



- > 1 for a consistent 



with the two-parameter conformal factor (5.19). Eq. (5.26) implies 

interpretation of the holographic RG flow. Panel (b) plots ^ 2L^^v " ^ function of 
the three-parameter conformal factor (5.27) with a = 5 x 10~^. 



^ for 



Of course, if one chooses a = 1, this profile reduces to the previous one in eq. (5.19).^^ 
Hence the AdS scales in the UV and IR limits are again given by eq. (5.20) with the 
new expression. However, in this case, the sharpness of the transition between the 
asymptotic UV and IR geometries is effectively controlled by both R and a. Given this 
new profile, we can readily extend the previous analysis to find in which parameter 
regime (i?, a) the entanglement entropy undergoes a phase transition. We do not 
present any details but a qualitative observation is that R^. becomes larger with smaller 
values of a. Now considering the FG expansion in this case, the metric function in 
eq. (4.8) becomes 

/(r) = (1 + 2ae~^^'^ + e"^^/^)""/' = 1 -^ae"^^'^ + .... (5.28) 

Hence we still have a = Luy / R and the lower bound in eq. (5.26) remains unchanged. 
Hence, as shown in figure 7b, we find that Rc can now satisfy this bound for > 3. 
Our expectation is that by further embellishing the form of the holographic RG flow 
proflle, we can also flnd realistic geometries, i.e., geometries satisfying eq. (5.26), which 
produce a phase transition in the entanglement entropy for = 2 as well. 



6. Holographic flow of c-function with GB gravity 

In this section, we return to examining the general behaviour of the c-function (3.4) in 
holographic RG flows but now where the bulk geometries are solutions of GB gravity 

^^With the choice a = 0, eq. (5.27) also reduces to eq. (5.19) upon substituting i? — > 2i? and 7 — > 27. 
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(2.9).^^ As noted in section 3.1, the calculation of holographic entanglement entropy in 
GB gravity requires that we extremize the entropy functional given in eq. (3.6) [49, 23]. 
Evaluating this functional for the strip geometry yields the expression in eq. (3.7). The 
extremal surfaces are again characterized by a conserved quantity which now takes the 
form given in eq. (3.8). The latter is most simply evaluated by considering the turning 
point of the extremal surface where r = and r = so that eq. (3.8) yields 



-{d-l)A{r^) 



(6.1) 



precisely as was found before for Einstein gravity. 

Recall that with Einstein gravity, we found a simple relation between and dS/M 
- see eq. (A. 12). In the following, we will show that this same relation extends to GB 
gravity. To simplify the discussion, we denote the integrand in eq. (3.7) as C and then 
Kd can be expressed as 

^ ^ dC , , 

C-—r. (6.2) 



Kd{() ~ dr 

Now we vary the entanglement entropy functional (3.7) with respect to the width of 
the strip £ to find 
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(6.3) 



Note that there is an extra overall factor of 2 above to since we are only integrating over 
half of the bulk surface, i.e., from the turning point (a;,r) = (0,rm) to the boundary 
(x,r) = ((^ — e)/2,rc). Now surfaces extremizing eq. (3.7) will satisfy 
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(6.4) 



Further eq. (4.5) still applies in the present analysis and so allows us to express (1 — 
de/dtj in terms of derivatives of the profile r(x, £). With both of these expressions, we 
are able to simplify eq. (6.3) to take the form 
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(6.5) 



^''We refer the interested reader to appendix B for a brief discussion describing the exphcit con- 
struction of such solutions. 
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Here the boundary term at x = vanishes because r = there. Eq. (6.2) was then used 
to simphfy the remaining terms with Kd{i). We note that precisely the same expression 
as above appeared in the analysis of dS/di with Einstein gravity - see eq. (A. 6). 

The next step is to show that the ratio 4 1^ has the same simple boundary limit as 
found previously in eq. (A. 11) with Einstein gravity. The analogous procedure would 
call for solving eq. (3.8) to find f. However, in the present case, we would find a cubic 
equation in and which would in general have six distinct roots. The relevant root 
would be that which in the limit A — )■ is continuously connected to the solution 
(A. 10) found with Einstein gravity. While it is possible to carry out this procedure 
analytically, it is not very illuminating. Rather we note that we are interested in the 
behaviour near the asymptotic boundary where the geometry approaches AdS space 
and the conformal factor takes the form A{r) = r/L. Now it is sufficient to expand 
near the boundary where is very large and the leading contribution to r becomes 

r ~ e*-/^K, (I-2A/0,) . (6.6) 

This equation is easily solved to yield x{r) near the boundary. Of course, the integration 
constant is chosen to satisfy the boundary condition x(r — i- 00) — i/2. Using this 
asymptotic solution for the profile of the extremal surface, it is easy to confirm that 



1 dr 
r di 



= -L (6.7) 

x=£/2 ^ 



as desired. Note that this relation is independent of the GB coupling A. 
Substituting eq. (6.7) into eq. (6.5), we arrive at 

dSGB 27rif'^"2 1 

(6.{ 



-p 



which precisely matches the expression (A. 12) found with Einstein gravity. It seems 
that this result is quite general. The first key ingredient is, of course, the conserved 
quantity (6.2). The other necessary ingredient is that the asymptotic geometry ap- 
proaches AdS space, which seems sufficient to produce the simple expression in eq. (6.7). 
Hence we expect that the expression (6.8) should be general to all cases with these two 
basic features. 

Next, we turn to the fiow of c-function (3.4) as we change the minimum radius 
of the extremal surfaces, as considered for Einstein gravity in section 4. Given eq. (6.8), 
our starting point for the c-function (3.4) is precisely the same as in eq. (4.20), i.e., 
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Hence using eq. (6.1), we find 



dcd ^ 2 7r((i - 
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Following our previous analysis, we express i in terms of with 
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However, in the present case, it is not possible to use the explicit root from eq. (3.8) 
for r and perform the integral. Hence we define 
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and use integration by parts in eq. (6.11) to write 

2 " -Ji^y^^^^) + / dr—F{r,rr, 
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(6.13) 



In eq. (6.12), we have chosen the above hmits on integration because the integrand van- 
ishes at the asymptotic boundary r = oo. Further, the dependence on in eq. (6.12) 
comes from r(x, r^). We should remind the reader that in the following discussion, our 
independent parameters are the profile r{x,i) and r^. Now differentiating eq. (6.13), 
we find 
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The first term above arises since dF{r, r^) /dr = A'{r)/r but note that this term should 
be evaluated slightly away from r = r^ since f = there. This potential divergence 
will be canceled below by a contribution which is revealed in the second term below. 
Now substituting eqs. (6.13) and (6.14) into eq. (6.10), we find 

47r(d-l)/3,£'^-2 
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and 
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In writing these expressions, we have used df/drm = {df / dKd){dKd/drm) , as well as 
eq. (6.1). Now considering the local expression for in eq. (3.8), we calculate dr/dK^ 
keeping r fixed. Similarly, we can differentiate eq. (3.8) with respect to x, which yields 
an expression involving f. We find that these two quantities are related by the following: 

rA'(^r + {d-l)Kd^^ = r + AXL^ A' A" Q (6.18) 



where 
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r2 (1 - 2 AL2A'2) + e2^(l + 4 AL2 A'^) 
We now use eq. (6.18) to express eqs. (6.16) and (6.17) as 
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where we have used J drf/f^ = — J dx d/dx{l/f) = —1/f. Now inserting these ex- 
pression (6.20) into eq. (6.15) and then integrating by parts, we arrive to following 
result 
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(6.21) 



Here we used the GB equations of motion, i.e., eq. (B.12), to replace A" by various 
components of the matter stress tensor. Note that the final result matches that in 
eq. (4.28) when A = 0. However, with A 7^ 0, it is clear that the null energy condition 
alone (which ensures T^t — T^r < 0) is insufficient to enforce a definite sign for dcd/drm- 
Rather we must also be able to make a clear statement about the positivity of the last 
factor in the integral, i.e., 



6 AL2A'2) + 



l-4AL2A'2g _ f2( 
1 - 2AL2 A'2 ~ (1 - 2AL2 A'2)(r2 (1-2 XL'^A''^) + e2^(l + 4 AL2 A'"^)) ' 



(6.22) 
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As described above, one could use eq. (3.8) to express in terms of the conserved 
charge Kd and the conformal factor A, however, the resulting expression for eq. (6.22) 
is lengthy and unilluminating. In the limit of small A we observe that eq. (6.22) becomes 



1 - A\L' A"Q 
1 - 2AL2 A'2 



2XL'A" + ---. (6.23) 



Hence it is not clear what simplification one might expect in eq. (6.21). However, this 
result is still suggestive in that it is easy to see that the right hand side is positive 
as long as A < 0. Unfortunately, examining the full expression in eq. (6.22), we see 
that this simple condition does not quite guarantee that this factor is positive. Thus 
while we have an expression for dca/drm in GB gravity, we are not able to make a 
simple statement of the conditions that are necessary to ensure that the c-function 
flows monotonically along holographic RG flows. 



7. Discussion 



With eq. (3.4), we constructed a simple extension to higher dimensions of the c-function 
(3.1) considered in ref. [2] for two-dimensional quantum field theories. As described 
in section 3, while the entanglement entropy itself contains a UV divergence, this ex- 
pression (3.4) is finite and, at conformal fixed points, yields a central charge that 
characterizes the underlying conformal field theory, as had been noted previously in 
[23, 46]. In section 4, we examined the behaviour of this c-function in holographic RG 
flows in which the bulk theory was described by Einstein gravity. In particular, we were 
able to show that the flow of the c-function was monotonic as long as the matter fields 
driving the holographic RG flow satisfied the null energy condition. As reviewed in 
section 2, the latter condition was precisely the constraint that appears in the standard 
derivation of the holographic c-theorem [18, 19, 20]. 

We observe that if the bulk geometry is such that it 'slightly' violates the null- 
energy condition over a 'small' radial regime, the integral in eqs. (4.19) or (4.28) would 
remain positive and hence the flow of our c-function would still be monotonic. That is, 
we only need the null energy condition to be satisfied in some averaged sense. Hence 
the null-energy condition is a sufficient but not a necessary condition for the monotonic 
flow of the c-function (3.4). Thus there is less sensitivity to the bulk geometry is the 
present construction of a holographic c-theorem using the entanglement entropy than 
in the original discussions [18, 19]. It is intriguing that when expressed as an integral 
over the boundary direction x, eq. (4.28) weights the contributions of the bulk stress 
tensor more or less equally for each interval 6x in the strip. However, when the integral 
is expressed as an integral over the radial direction, the integrand includes an extra 
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factor or 1/r, which diverges at the minimum radius Vm of the holographic surface (but 
the integral remains finite). Hence in the holographic flows, the c- function responds 
sensitively to changes in the geometry near this radius in the bulk - a result that can be 
seen in the explicit flows discussed in section 5. Hence given the holographic connection 
between radius in the bulk and energy scales in the boundary theory, it seems clear 
that the flow of this c-function is most sensitive to the lowest energy modes probed by 
the entanglement entropy. 

Our result for the monotonic flow of c-function in section 4 refers to the derivative 
dcd/drjn, i-G-, changes in Cd as we change the minimum energy scale probed by the 
entanglement entropy. To describe the flow of q completely in terms of the boundary 
theory, we would actually like to establish dca/di < 0, i.e., the c-function decreases 
monotonically as we increase the width of the strip for which the entanglement entropy 
is evaluated. In this case, we would be using the width £ as a proxy for the relevant 
energy scale along the RG flow. The desired result can be established in the present 
holographic framework, however, as discussed in section 5, one must be careful to 
restrict attention to the physical saddle points in evaluating the entanglement entropy. 
We showed there that extremal surfaces can arise for which di/drm > and hence 
dcd/d£ > 0. However, these saddle points do not contribute when one evaluates the 
holographic entanglement entropy with eq. (3.5) since they are never the minimum area 
surface. Rather the appearance of these 'unstable' saddle points signals a first order 
'phase transition' in the entanglement entropy. As a result, Cd drops discontinuously at 
some critical value it of the width of the strip and the monotonic 'flow' of the c-function 
is preserved. 

While we have only illustrated this behaviour with specific examples in section 4, 
it seems clear that the physical entanglement entropy will never be determined by such 
saddle points. In particular, if we are studying a holographic RG flow between two 
AdS geometries, we will always find di/drm < when is well into either of these 
two asymptotic regions. Hence as argued in section 5.2, if extremal surfaces arise for 
which di/drm > 0, it indicates that there are a number of competing saddle points in 
the corresponding regime. First i{rm) is a single-valued function since the conserved 
charge (3.8) dictates that there is a unique surface for each value of r^. Hence if 
we assume this is a smooth function, there will always be (at least) three competing 
saddle points when di/dr^, > 0. It then becomes inevitable that there will be a phase 
transition in the corresponding regime of i. Further we note that 

dS _ H^Cddi_ 

d^ ~ l^Jdd^r 

The first two factors above are positive definite and hence the sign of dS/dr^ is con- 
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(7.1) 



trolled entirely by d(i/drm- Given this result, it is straightforward to argue that the 
behaviour illustrated in figure 4 is in fact generic. That is, the phase transition goes 
between the two branches for which M/drm < 0. Hence we have argued that given 
dcd/drm > 0, it also follows that dcd/dl < for RG fiows dual to Einstein gravity. 

One may be concerned that the phase transitions noted above are an artifact of 
choosing a background geometry in the bulk which is unphysical in some way. However, 
with our analysis in section 5 and appendix C, we argued that the phase transitions 
can arise for holographic backgrounds that have a natural interpretation as an RG fiow 
in the boundary theory, but also for backgrounds where the interpretation seems to be 
more exotic. While this interpretation was explicitly shown to apply in examples of 
phase transitions with the boundary dimension c? > 3, constructing further examples 
to extend this result to = 2 does not seem difficult. However, we note that these 
phase transitions are undoubtedly effect of the large limit which is implicit in our 
constructions. However, it may still be that similar behaviour, i.e., rapid transitions 
in the entanglement entropy, persists in the RG flows of more conventional physical 
systems. In any event, it would be interesting to better understand these phase tran- 
sitions in the holographic systems. Such a transition seems to indicate that quantum 
correlations in underlying degrees of freedom change dramatically at some particular 
energy scale in the RG flow. 

It is worthwhile noting that phase transitions in the holographic entanglement en- 
tropy of the kind found here and in [26] for RG flows also arise in a variety of other 
holographic constructions. The simplest example is to consider the case where the 
entangling surface contains two disjoint regions. When the two regions are relatively 
close together, saddle point determining the holographic entanglement entropy will be 
a single connected bulk surface. However, as the two regions are moved apart, there is 
a phase transition to a second saddle point consisting of two separate bulk surfaces [47] . 
A similar phase transition was also found in considering the holographic entanglement 
entropy of the strip geometry for a bulk background corresponding to a confining phase 
of the boundary theory [54, 55]. There is a strong similarity between the results for 
these confining theories and the present RG fiows since the phase transition again arises 
as the width to the strip passes through some critical value and results in a discontin- 
uous drop in the central charge q. Further, in both cases, the phase transition can be 
interpreted as being produced by a rapid and drastic restructuring in the correlations 
of the low energy degrees of freedom (in comparison to high energy correlations). Sim- 
ilar results were also found for other entangling geometries, i.e., a circular surface in 
three-dimensional confining boundary theory [56]. Finally similar phase transitions in 
the entanglement entropy have also been found in holographic superconductors as the 
temperature is varied [57] and in the time evolution of holographic quantum quenches 
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[58]. 

In section 6, we considered extending our results to holographic models where the 
gravitational theory in the bulk is Gauss-Bonnet gravity (2.9). While it is straightfor- 
ward to construct an expression (6.21) for dcd/drm in GB gravity, it is evident that the 
null energy condition is not sufficient to guarantee a monotonic flow of the c-function. 
Unfortunately, eq. (6.21) does not lend itself to a simple statement of the conditions that 
would necessary to ensure that the c-function flows monotonically along holographic 
RG flows in these models. Further insight into this question may be provided by ex- 
amining explicit holographic RG flows. In section 5, we assumed that the holographic 
backgrounds were solutions of Einstein gravity and hence the entanglement entropy is 
determined by eq. (3.5). We could just as easily assume that the same backgrounds 
are solutions of GB gravity and examine the behaviour of the c-function defined by 
eq. (3.6). In particular, it would be interesting to see if there are violations of the 
monotonic flow of the c-function in certain parameter regimes. 

Of course, it may not be a surprise that the monotonic flow of the c-function (3.4) 
is not directly connected to the null energy condition in GB gravity. As described in 
section 2, an important feature of this theory is that at conformal flxed points, the 
dual boundary theory has two distinct central charges, given in eqs. (2.13) and (2.14). 
Using the null energy condition, ref. [20] established that the charge denoted a*^ would 
satisfy a c-theorem in these holographic models. However, in section 3, we found that 
the c-function (3.4) actually corresponds to a nonlinear combination of both central 
charges. Hence as we noted at the outset, it was improbable that a simple holographic 
c-theorem could be established for GB gravity with the present construction. Setting 
holography aside, it is known that for four-dimensional quantum fleld theories, there 
is no possible (linear) combination of the two central charges, c = Ct and a = ajj, that 
can satisfy a c-theorem other than a alone [8]. 

Of course, GB gravity only provides an interesting extension of the usual holo- 
graphic framework for d > A. For smaller values of d, the curvature-squared interaction 
(2.10) does not contribute to the gravitational equations of motion because of the topo- 
logical origin of this term. It may be of interest to study the behaviour of our c-function 
for other gravity theories with higher curvature interactions for d = 2 and 3. Inter- 
esting families of holographic models were considered with higher curvature theories 
of the three-dimensional gravity in [53]. A deflning feature of these theories was that 
the dual d = 2 boundary theory should exhibit a c-theorem. Hence these models may 
provide an interesting holographic framework to examine the RG flow of C2- However, 
the work of [2] indicates that this flow must be monotonic for any unitary and Lorentz 
invariant quantum fleld theory and so conflrming this result here would really be a test 
that these holographic models deflne reasonable boundary theories. 
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Our construction of the c-function (3.4) can be applied quite generally, i.e., outside 
of the context of holography. While the RG flow of Cd is not expected to be monotonic 
in a generic setting, we observe that the flow can be constrained somewhat following 
the analysis of [2]. In particular, let us define a new (dimensionful) function of the 
following form in arbitrary d: 



Of course, for d = 2, we have C2 = C2. In any event, we can apply directly the analysis 
of [2] to show that dcd/M < for any d. There, the authors considered unitary and 
Lorentz invariant field theories and used sub-additivity of the entanglement entropy 
to show that C2 is a monotonically decreasing function as £ increases. In particular, 
they considered two specific surfaces h and c with a relative boost, as shown in figure 
8. Further a and d chosen as constant time surfaces in some frame so that they are 
Cauchy surfaces whose causal development corresponds to the intersection and union 
of the causal development of the original two boosted surfaces. By construction the 
surfaces a, h and c just touch the boundary of the causal development of d on either 
end. Two important observations are: First, if we are evaluating the entanglement 
entropy for these segments in the Lorentz invariant vacuum state, then it should only 
depend on the proper length of the corresponding interval. Second, the entanglement 
entropy of any of these surfaces will be the same as for any other Cauchy surface of 
the corresponding domains since the time evolution is assumed to be unitary. Now the 
authors of [2] show that sub-additivity of the entanglement entropy of these regions 
imposes the following relation 



where they introduce the ratio of the proper lengths a = c/a = d/b > 1. Now if the 
relative boost is taken to be small, this relation implies that £dS/d£ = C2/S = C2/S must 
decrease with increasing i. Now these elegant arguments can also be applied without 
change in considering the entanglement entropy for the strip geometry in arbitrary 
dimensions, as long as the boundaries of the strip are orthogonal to the boost direction. 
That is, the segments a, b, c and d in figure 8 now represent the orthogonal cross-section 
of four specific strips. Hence with the previous arguments, we will recover the relation 
(7.3) for arbitrary dimensions. This implies that with the strip geometry in arbitrary 
dimensions, we have dcd/di < — the same observation was made in [54]. While 
we cannot conclude that the original c-function (3.4) must decrease with increasing i, 
using eq. (7.2), we find an upper bound on the rate at which q could increase: 




Pd ndS 



(7.2) 



H'^-^ di ■ 




(7.3) 



dcd 
d£ 



< {d-2)cd. 



(7.4) 
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This bound applies to any c?- dimensional quantum field theory, subject to the provisos 
of Lorentz invariance, unitarity and sub-additivity of entanglement entropy. 




Figure 8: Here the boosted surfaces b and c are drawn with time and space on the vertical 
and horizontal axes, respectively. The casual development of the strips d and a contain the 
union and intersection of the causal development of b and c. The causal development for each 
of these surfaces is outlined with dashed lines. 

An important feature of the entanglement entropy (3.2) for the strip geometry is 
that it contains a single divergent term which is independent of the width i. As a 
result of this simple structure, the single derivative in eq. (3.4) produces a UV finite 
or regulator independent c-function. Now for a general smooth entangling surface, the 
structure of the divergent contributions to the entanglement entropy at a conformal 
fixed point is more complicated, i.e., 

^ ^ Pd^ Pd~4 (f+po + 0{6) for odd rf, 

Sd-^^6''~^ \^+Polog{l)+0{l) for even d, ^ ' 

where 6 is the short distance cut-off (and i is some IR scale). Of course, to produce a 
dimensionless entanglement entropy, the coefficients pn must have dimension length". 
At a conformal fixed point, this dimensionful character is provided by scales arising from 
the geometry of the entangling surface and the background spacetime. For example, 
the first term yields the celebrated area law with Pd~2 oc A, where A is the area of 
the entangling surface. If we move away from a conformal fixed point, there will be 
additional dimensionful couplings in underlying theory and the divergence structure of 
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the entanglement entropy can become substantially more complicated. For example, in 
a holographic setting, the additional complications of relevant operators were illustrated 
in [24]. 

Of course, the conjecture by [20] proposes that a central charge, which is identified 
using the entanglement entropy, should satisfy a c-theorem (1.1) in higher dimensions 
. In particular, the prescription there specifies that the entangling surface should be 
a ((i-2)-sphere of radius R in flat space. However the universal central charge was 
identified with the dimensionless coefficient po which only appears at subleading order 
in eq. (7.5) for d > 3. Further since R is the only scale in the problem, one expects all 
of the preceding coefficients are nonvanishing with p„ = p„ i?" for some dimensionless 
coefficients p„. Hence a single derivative with respect to the radius of the sphere 
will not isolate the desired coefficient in the entanglement entropy. Certainly a more 
sophisticated construction would be needed to remove all of the potentially divergent 
terms along a general RG flow. In fact, precisely such a construction was recently 
proposed in [26]. It would be interesting to investigate this new proposal along the 
lines of the analysis in sections 4 and 6. In particular, it would encouraging if the 
monotonic flow of the c-function identifled with this construction or a similar variant 
could be directly related to null energy condition in holographic RG flows. 
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A. The c-function with Einstein gravity for d > 3 

In this appendix, we describe the details of the derivation of eq. (4.20), which is our 
starting point for the discussion in section 4 of the holographic flow of the c-function 
in higher dimensions. The analysis here is largely an extension of that given there for 
c? = 2 to an arbitrary d. For a d-dimensional boundary fleld theory, the holographic 

^^The applicability of this geometry to the c-theorem conjectured in [20] follows from the results of 
[21]. 
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entanglement entropy of the strip is given by 



S 



dx e('^-2)A(r) ^-2 ^ g2A(r) 



and the conserved charge (3.8) reduces to 



(A.l) 



(A.2) 



Further, for a particular extremal surface, = e~^^~^^^^^"^^ since at the turning point 
we have r(0, i) = and f(0, i) = 0. 

Now we differentiate the entropy (A.l) with respect to the width i while holding 
the radial cut-off Vc fixed. In doing so, we keep in mind that both the radial profile for 
the surface and the cut-off in x are implicitly functions of the width i, i.e., r = r{x,i) 
and e = e{i). Hence we find 



_ 4 7rij-'^-2e('^-2)A(r) 



de 



dr 

m 



(A.3) 



x=(i-e)/2 



where f = dr{x,i)/dx. In simplifying dS/di to produce the expression above, we have 
used r(0, £) = to eliminate the boundary terms at x = 0, as well as removing the bulk 
variation using the equation of motion which follows from eq. (A.l), 



r-A' {dr' + {d-l)e^^) 



0. 



(A.4) 



The cut-off e is defined by the relation: r(x = {i — e)/2,£) = re- We can vary i in the 
latter expression while holding rc fixed, as in the variation in eq. (A.3), to find 



r[x, ■ 



de 
dl 



dr{x, tj 
di 







(A.5) 



x=ie-t)/2 



This result was already presented in eq. (4.5). Using this expression and also eq. (A.2), 
we can now simplify eq. (A.3) to take the form 



dS 



inH'^-^ 1 dr 
r di 



(A.6) 



x=(£-e)/2 



Further progress requires that we consider the bulk geometry in the far UV region, 
which we assume approaches AdS space asymptotically with 



A{r) = r/L. 



(A.7) 
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Now recall that eq. (A. 2) yields a simple expression for dx/dr, as presented in eq. (5.2). 
With the simple conformal factor (A. 7), this equation is easily integrated to yield 



X 



L e-*/-^ 
dKd{l) 



d 



3d -2 e-^C'i-iWi 



2'2(rf-l)'2(c/-l)' Kdiiy 



(A.8) 



where we are imposing the boundary condition that x — i/2 asr — t- oo. We note 
that the same hypergeometric function appears in eq. (4.25). This is not a coincidence 
because in the present case A' is simply a constant and so we are essentially performing 
the same integration here as in eq. (4.25). Further note that this complicated expression 
reduces to eq. (5.3) for d = 2. Taking the partial derivative of eq. (A.8) with respect 
to we find 



dr{x, 



di 



2 

dKd 

{d-l)Kd di 



^l)r/L 



(A.9) 



L 



- e-2(d-l)r/L 



dK, 



d 



2' 2(rf- 1)' 2(rf- r 



U-2 e-2{'^-i>'/^ 



KAi)^ 



Now using eq. (A. 2) with the conformal factor (A. 7), we find 



dr/L _ g-2(d-l)r/L 



(A.IO) 

Then we consider the ratio of these two expressions and take the limit x ^ i/2 and 
r — > oo to find 

Essentially this result indicates that with the limit — oo, de/di vanishes in eq. (A. 5). 
Now substituting eq. (A. 11) into eq. (A. 6) yields 



dS 



2-kH' 



d-2 



We use this simple form of derivative of entanglement entropy to produce the expression 
for the c-function in eq. (4.20). 



B. RG flow solutions for GB gravity 

We are interested in holographic RG flow solutions for the action (2.9) in which GB 
gravity is coupled to a scalar field, i.e., 



2ii~^ 



'' \d-2)id-,f ^-V'^^'-^^^\ 



(B.l) 
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where is given by eq. (2.10). First, we will examine the equations of motion for holo- 
graphic RG flows in some detail. This will allow us to explicitly verify the result used 
in eq. (2.18) to prove the holographic c-theorem for GB gravity [20]. Then we discuss 
a simple approach to explicitly construct analytic solutions for GB gravity describing 
holographic RG flows. This construction is a simple extension of the 'superpotential' 
approach developed for Einstein gravity in [59, 60]. This approach was already ex- 
tended to GB gravity with four boundary dimensions in [61] and here we provide the 
generalization to arbitrary d > 4. 

In the action (B.l) above, the cosmological constant term has been absorbed into 
the scalar potential and we assume that V{(j)) has various critical points where the 
potential energy is negative as in eq. (2.2). As described in section 2, for each of 
these critical points, there is an AdS vacuum solution where the curvature scale is 
given by = L^/foo where L is some canonical scale appearing in the potential and 
the curvature squared interaction while the (dimensionless) constant f^o is given by 
eq. (2.12). To consider solutions describing holographic RG flows, we begin by writing 
the scalar and gravitational equations of motion as 

V20-^ = O, (B.2) 

1 AL^ 
Rab - ^Rgab + (rf_2)(rf_3)^'^^ = ^'^^ ' ^^-^^ 

where ^ 

Hab = RacdeRb'^'^^ — '^RacRb'^ — '^RacbdR'^'^ + RRab — -^^iQab ■ (B-4) 

Further the stress tensor for the scalar field is given by 

1„ , „ , 1 /I ,,2 



Tab = ^da<Pdt<P--gab\^-{d<pr + V{<P)j . (B.5) 

As in the main text, we consider the following ansatz for the metric: 

ds^ = e2^(") r]^j dx'dx^ + cir' , (B.6) 

For nontrivial RG flows, we also include a simple ansatz for the scalar: = 0(r). In 
particular then, this ansatz maintains Lorentz invariance in the boundary directions. 
Now with these metric and scalar ansatze, there are two nontrivial components of the 
gravitational equations (B.3): 

d(d-l)(A'2-ALM'^) = 2T/, (B.7) 
2(rf-l)(l-2ALM'2)A" + rf(rf-l)(A'2-ALM'^) = 2Ti*. (B.8) 
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Again 'prime' denotes a derivative with respect to r. Using eq. (B.5), we find the 
following components of stress tensor 

T/ = \{<Pr-\v{<P). (B.9) 

T,^ = -6,^(^-{<P'f + \v{<P)^ • (B.IO) 

With the present ansatz, the equation of motion for the bulk scalar (B.2) becomes 

+ ^ = 0. (B.ll) 

00 

Of course, the three equations of motion above are not all independent. For example, 
one can derive eq. (B.8) by differentiating eq. (B.7) and then substituting in eqs. (B.7) 
and (B.ll). 

Taking the difference of eqs. (B.7) and (B.8), we find 

Tt' - T/ = (rf - 1)(1 - 2 A'2)A" , (B.12) 

which is the result used to establish the holographic c-theorem for GB gravity in 
eq. (2.18). Eqs. (B.7), (B.8) and (B.12) are written for GB gravity coupled to a gen- 
eral matter field Lagrangian and as long as the matter sector satisfies the null energy 
condition, the combination of components of Tab in eq. (B.12) are negative (or zero). 
In the particular case considered here, i.e., the action (B.l) with a scalar field, we find 
T^t — T^r = ~(0')^/2 and so the sign of eq. (B.12) is obvious. 

Given eqs. (B.7-B.11), a simple set of explicit solutions can be constructed by 
extending an approach developed for Einstein gravity in [59, 60]. The key idea is to 
consider a special class of scalar potentials that can be defined in terms of a 'superpo- 
tential' and then express the solution in terms of this superpotential. This construction 
was extended from Einstein gravity to GB gravity in five dimensions {i.e., d = A) in 
[61] and here we provide the generalization to arbitrary d > 4. First, we write the 
scalar potential in terms of a superpotential W{(j)) as follows: 



V{^) = 2{d 



1)2 (^^y {l-2XL^W'Y -d{d-l)W'il-XL'W^). (B.13) 



With a potential of this form, the equations of motion above can be re-expressed as 
first order equations: 

0' = _2(rf-l)(l-2AL2W^2)^, 

00 

A' = W. (B.14) 
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Given these first order equations, we may now solve for tlie metric (B.6) and tlie scalar 
profile in quadratures. We note that the same equations appear in the construction of 
domain wall solutions for 'new massive gravity' in higher dimensions [62].^^ The action 
studied there differs from the GB action (B.l) by the addition of an action proportional 
to the square of the Weyl tensor, i.e., — (■^_2)^rf_3) CabcdC"'^^'^. This additional interaction 
does not influence the equations of motion in the present setting because the metric 
ansatz (B.6) for the holographic RG flows is conformally fiat. 



C. Scalar potentials for section 5.2 



In section 5.2, we analyzed the flow of the entanglement entropy in various bulk geome- 
tries which were defined by specifying an explicit conformal factor, e.g., as in eq. (5.19). 
Here we would like to show that these profiles can arise as solutions of Einstein grav- 
ity coupled to a scalar field with an appropriate potential. In particular, we use the 
'superpotential' approach described at the end of the last appendix, of course, after 
setting A = 0. With the latter simplification, eqs. (B.13) and (B.14) reduce to 



V{^) = 2{d-lY 
(f)' = -2{d- 1) 



6W 

6W 

1^' 



d{d-l)W\ 
A' = W. 



(C.l) 
(C.2) 



Now using these equations, we would like to construct V{(f)) given A{r) 
First from eq. (C.2), we find 



A" 



6W 



64> 



-2{d 



6W 



k'\2 



2{d 



Hence if we let 
we can write 



at the critical point in the UV, i.e., (pir — )• oo) 



dr(j)'{r) = 



dr 



2{d-l)A"{f) 



1/2 



(C.3) 
then 

(C.4) 



Similarly combining eq. (C.l) with eqs. (C.2) and (C.3), we may write the value of 
potential along the flow as 



V{r) 



id-l)A' -d{d-l){A 



l\2 



(C.5) 
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Given these two expressions, one can easily make a parametric plot of the potential V 
as a function of (p, at least over the range of the scalar covered in the holographic flow. 
Alternatively, if eq. (C.4) can be integrated analytically and inverted, i.e., one can 
write r = r(0), then eq. (C.5) will yield an analytic expression for V{(j)). We use both 
of these approaches to describe the potential corresponding to the conformal factors 
presented in section 5.2. 

Let us begin with the conformal factor given in eq. (5.19). In this case, it is 
straightforward to integrate eq. (C.4) with the result being 



tan I — 

\2 0UV 



with = TT ^/2^{d - 1) . 



(C.6) 



We also note that at the far IR of the holographic flow, — j- as r — ?■ — oo. Given this 
expression (C.6) and the conformal factor (5.19), it is straightforward to calculate the 
superpotential W{(f)) using eq. (C.2), 



1 7 
W{(f>) = 2^ R^'^^^ 



(C.7) 



and the potential V{(f)) using eq. (C.5), 

m = ^sm^^^V^(^ 



i?2 



1 

L 



1 
R 



cos vr 



(C.^ 



Hence we have produced a analytic result for the scalar potential necessary to produce 
the bulk geometry with A{r) as in eq. (5.19) as a solution of Einstein gravity. In figure 
9, the curve with a = 1 illustrates the behaviour of the potential (C.8). Figure 10 also 
shows the potential for various values of the parameters. 

Now let us consider the conformal factor given in eq. (5.27). Using the same 
equations as before, we now find 



2^-f{d-l) V 1 + ^ cosh(2r//2) 



dr 
V{r) 



R 

27(rf-l) 



a + cosh(2r/i?) 
- a cosh(2r/i?) 



(C.9) 



{a + cosh{2r/R)y 



d{d 



1 7 sinh(2r/i?) 
Z ^ ^ (T + cosh (2r/i?) 



In this case, we were not able to analytically integrate 0'. However, it is straightforward 
to perform a numerical integration and produce a parametric plot of the potential V as 
a function of 0, as shown in figure 9. In this figure, the boundary dimension was chosen 
to be (i = 4 and the potential parameters are fixed as L = .40, R = .55 and 7 = .825, 
while a is varied as indicated on the plot. The plot shows that qualitatively this 
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Figure 9: A parametric plot of y(0)/|\^v| versus (f) for eq. (C.9) with various values of 
parameter a. The other parameters in the potential are set to L = .40, R = .55 and 7 = .825, 
as well as setting d = 4. 

potential has the same shape as the analytic potential found for the previous conformal 
factor in eq. (C.8) - recall that in fact, this new potential will agree with eq. (C.8) when 
cr = 1. However, the new potential is becoming steeper as the parameter a becomes 
smaller. The plots only show the potential over the range relevant for the holographic 
RG flow, i.e., from = = to 0uv However, the full potential would be symmetric 
about both and 0uv and so it would be periodic with a period A0 = 2(0uv — (pm)- 
Further we note that for the parameters chosen in figure 9, the entanglement entropy 
undergoes a phase transition for a = .0005 but there is no phase transition for either 
a = .2 or 1. 

While we were not able to produce an analytic expression for the full potential 
from eq. (C.9), it is still possible to consider a perturbative expansion around the UV 
critical point, i.e., 



where 6(f) = </) — </)uv One can easily verify that there are no odd powers of 6(j) appearing 
in this expansion. As an example, we note that the mass parameter above would be 
calculated as 



r (50) = Vuv + 2 ^uv ^0' + - Auv 50' + ■ ■ ■ , 



(C.IO) 





The results for the first few parameter is the expansion (C.IO) are: 

d(d — 1) n 1 d 

Kiv = -r^ - 



t2 ' UV P2 r p 

. ^d LuvR + d{2 - 7a^) LjrR - 8(1 - 2a^) LirLuv , . 

2{d-l)aHLuv-L,n)R^ ' ^^-'"^ 

where we are using the expressions in eq. (5.20) for Luv and Ljr. Of course, the value 
of Vuv above corresponds to the expected value of the cosmological constant in the 
asymptotic AdS geometry. We note that the mass parameter is independent of 
the extra parameter a appearing in the conformal factor (5.27). Hence these first two 
parameters precisely match those in the analogous expansion of the analytic potential 
in eq. (C.8), as is readily verified. As shown in eq. (C.12), the additional parameter a 
first makes its appearance the quartic coupling Auv In this case, one can still verify that 
Auv(cr = 1) matches the analogous quartic coupling found in expanding the analytic 
potential (C.8) about (p = 0uv To provide some qualitative insight for this quartic 
coupling, we add that typically {e.g., for the parameters chosen in figure 9), Auv is 
negative for cr ~ 1 but it becomes positive for small values of cr (and diverges as 
a 0). 

At this point, we would like to consider the holographic interpretation of these 
geometries in more detail. Recall that the standard description begins with a discus- 
sion in the boundary theory where a UV critical point is perturbed by some relevant 
operator O and the latter triggers an RG flow to a new critical point in the IR. A 
natural assumption in this discussion is that O is relevant, i.e., that it has a conformal 
dimension A < d. Now in the gravity description, this operator is dual to the bulk 
scalar and the conformal dimension is related to the scalar mass by the standard 
formula m^^L^y = A(A — d) [63]. Inverting the latter relation yields two roots for A, 
i.e., 

, jl 

2 V 4 

where the standard choice corresponds to A = A+. The scalar will have two indepen- 
dent solutions asymptotically [63], 



A± = -±yl-+ml^Lly, (C.13) 



where the coefficient of the more slowly decaying solution corresponds to the cou- 
pling for the dual operator while 0*"^' is proportional to the expectation value (O). 

Now turning to the scalar mass in eq. (C.12), we observe that the result can be 
written as m^^Lfjy = a{a — d) where a = Luy/R. Further by combining eq. (5.28) 
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Figure 10: Plot of y((/))/|Kjv| 

versus (/>/(/>uv for eq. (C.8) with various values of parameter 
LjjyjR. The parameters in the potential are chosen to fix L\jy = 1 and L//j = .25 and we 
set d = 4. 

and (j)' oc \J —N' from eq. (C.3), we see that a is precisely the exponent controUing the 
asymptotic decay of 50. Now there are four regimes of a which we consider separately: 
(a) < a < c//2, (b) c//2 < a < (c/ + 2)/2, (c) (rf + 2)/2 < a < and (c) d < a. In 
the interval (a), the holographic interpretation of the flow geometry precisely matches 
that described above. That is, a = A__ with A_|„ = d — A_ < d. Hence the dual 
operator is relevant and leading contribution in the asymptotic decay of 6(j) reflects 
the fact that corresponding coupling is nonvanishing in the boundary theory. In the 
intervals (b) and (c), we have instead a = A4. < d and hence the interpretation is that 
O is relevant but the boundary coupling for this operator vanishes. Rather holographic 
interpretation of the bulk solution is that O has a nonvanishing expectation value in the 
UV, which then triggers the RG flow to an new IR fixed point. Of course, this appears 
to be a somewhat unconventional description of RG flows. We have distinguished 
the interval (b) with | < a < | + 1 because in the regime, one can imagine the 
standard interpretation of the bulk solution still holds for the 'alternate quantization' 

^"Similar flows arising in a supergravity construction were studied in [64]. However, the flows 
described there would all fall into the interval (b). 
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of the holographic theory [65]. In the latter approach, the roles of 0*"^ and 0^+-' are 
reversed and the dimension of the boundary operator is given by A = A_. Note 
that in these first three intervals (a), (b) and (c), we have m^y < and m^y reaches 
its minimum value at a = d/2, where it coincides precisely with the Breitenlohner- 
Freedman bound [66], i.e., m^^Lljy = —d'^/4. Finally in the interval (d), we have 
a = > d and hence in this regime, O is no longer a relevant operator. Rather 
we have an irrelevant (or marginal) operator again with a nonvanishing expectation 
value which triggers the flow to a new fixed point in the IR. Note that with a > d, 
rrtyy > and the scalar potential must have a third extremum 03 between the UV 
and IR critical points where ^(03) > V^uv^^ This extra extremum is distinguished 
from the UV and IR critical points since by construction, the latter are critical points 
of the superpotential, whereas the new critical point satisfies 51^/501 ^=,^3 = but 
(5W^/(50|^=^3 7^ 0. The changing structure of the potential as Ljjv / R varies through 
these different regimes is illustrated in figure 10. As already noted, the interpretation 
of the RG flows in cases (c) and (d) in the boundary theory seems somewhat exotic. 
Hence one must worry that the underlying holographic model for these constructions 
is unphysical in some way. For example, it could be that for a consistent boundary 
CFT, once the dimension A of the operator is fixed, the quartic coupling in eq. (C.12) 
must be constrained in some way, along the lines of various bounds found in [34] or 
[67]. Alternatively, it could be that these background solutions are simply unstable in 
the corresponding parameter regime. In any event, in the main text, we focus our 
attention on the models in the intervals (a) and (b) for which the holographic RG flows 
have a conventional interpretation - see the discussion around eq. (5.26). 

To close this section, we note that one can also consider a perturbative expansion 
of the scalar potential as in eq. (C.IO) but about the IR critical point, i.e., 

V(M) = Kr, + ^ m2^50' + ^ A:h,50' + ■ • • , (C.15) 

where (50 = — 0ir with 0ir = 0(r — —00). As before, there are no odd powers of 50 
appearing in eq. (C.15). The first few parameters in this IR expansion are: 

d{d — 1) 2 1 



. -d{2- 7a^) LuvR - 3d LirR - 8(1 - 2a^) LirLuv 

2{d-l)aHLuv-Ljn)R^ ' ^^'^^^ 

As in the UV expansion, we see that the parameter a first appears here in the quartic 
coupling Air. Hence one can easily verify that V[r and m'f^ precisely match the corre- 

^^In fact, for Ljjv / R sufficiently large, one finds V{(f>-s) > 0. 
^^We thank Alex Buchel for discussions of these issues. 
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spending couplings in the IR expansion of the analytic potential in eq. (C.8), while Ajr 
matches the quartic coupling in this expansion when a = 1. Note that m'^^ > and 
so the dual operator is always irrelevant at the IR critical point. The behaviour of the 
quartic coupling is qualitatively the same as described above for Air, i.e., it is typically 
negative for a ~ 1 and positive for cr ~ 0. 
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